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I  INTRODUCTION 


A  GENERAL 


The  first  radio  reflections  from  the  moon,  using  conventional  radar 
techniques,  were  recorded  at  a  frequency  of  111  5  megacycles  at  moonnse  ar.d 
moonset  by  the  U  S  Army  Signal  Certs  ^  A  more  complete  description  of  the 
Signal  Corps  results  ar~  ted  by  je  Witt  arid  Stoic  la  *  In  order  to  obtain  the 

required  sensitivity,  a  receiver  bandwidth  of  about  5C  ep*s,  transmitter  pulse 
length  of  0-25  second  and  a  transmitter  peak  power  cf  three  kilowatts  were  used. 
The  experimental  results  shoved  tnat  the  returned  echo  amplitude  was  often  times 
less  than  the  theoretic  illy  computed  ere  and,  on  occasions,  could  not  be 
detected  In  addition,  the  signal  amplitude  r  ’  '  '°nt  .  ?/l-‘  :rr-u  :  inctuV  ions 

having  period!  cf  several  minutes  which  w^re  attributed  to  anomalous  ionospheric 
refraction 


At  approximately  the  same  time  j  the  Signal  Corps  observations,  Bay 
also  reporfe  that  radar  contact,  had  been  made  with  the  moon  at  a  frequency  of 
120  megacycle  These  measures  nts  indicated  a  pev^r  reflection  coefficient,  for 
the  noon’s  surface,  of  the  order"  cf  0  1 

Kerr,  Sham  and  Higgins  cf  Australia  performed  th°  moon  reflection 
experiment  at  21  5  megacy  V~  during  moo nr i se  and  moonset,  mainly  to  study  the 
characteristics  of  low  ingle  propagation  through  the  ionosphere  Their  records^ 
indicated  a  rather  slow  amplitude  fluctuation  rate  which  also  was  attributed  to 
ionospheric  refraction  phenomena.  The  rapid  amplitude  fluctuations  having 
periods  of  a  few  seconds  were  explained  by  assuming  that  the  moon  s  surface 
consisted  of  many  random  sea  ter^rs  ir.  relative  motion,  brought  about  by  the  moon 
libration  This  phenomenon  tcgs  then  with  rh^  evidence  that  the  received  pulse 
was  stretched  out  several  ♦ im^s  vh'  length  of  the  (one  millisecond  transmitted 
pulse,  showed  that  the  meet  was  a  "rough"  r^fl^etcr  at  their  frequency 


1 


The  first  successful  transmissions  at  UHF  using  the  moon  as  a  relay 
station  were  accomplished  by  Sulzcr,  Montgomery  and  Gerkt  CW  radio  signals  at 
a  frequency  of  4l8  megacycles  were  relayed  from  Cedar  Rapids,  Iowa  to  Sterling, 
Virginia  It  was  found  that  the  signal  was  subject  to  severe  fading,  its 
amplitude  varying  from  the  receiver  noise  level  to  occasional  peaks  as  high  as 
10  db  aoove  the  noise. 

7 

Measurements  of  radar-lunar  echoes  by  Evans  at  120  megacycles,  by 

8  9 

Trexler"'  at  198  megacycles  and  by  Yaplee,  Bruton,  Craig  and  Roman'  at  2860 
megacycles  revealed  that  radar  reflection  takes  place  within  a  small  area  at  the 
center  of  the  moon's  surface,  the  radius  of  which  is  approximately  one-third  the 
radius  of  the  moon  These  results  confirmed  that  the  moon  appeared  to  be  a 
quasi-smooth  reflector  at  radio  wave  frequencies. 

The  characteristics  of  lunar  echoes  such  as  libration  fading  of  the  order 
of  two  to  three  cps  and  pulse  lengthening  to  approximately  one  lunar  radius  have 
recently  been  noted  by  investigators  operating  in  the  k()0  to  UUO-rnc  frequency 
range. ^9  ^  The  latter  phenomenon  indicated  the  possibility  that  the  rear 
portion  of  the  moon  behaved  as  a  rough  scatterer 

The  employment  of  radar  techniques  co  reflect  signals  off  the  moon 1 s 

surface  for  studing  the  characteristics  of  the  ionosphere  by  means  of  the 

12  l?  lk 

Faraday  effect  was  first  reported  by  Browne  and  Evans.  9  Their  equipment 

consisted  of  a  10-kv  peak  power  transmitter,  operating  with  a  30-millisecond 

pulse  length  and  a  pulse  repetition  frequency  of  0  6  pulses  per  second,  and  an 

antenna  whose  size  was  approximately  250  square  meters.  The  degree  of  rapid 

echo  amplitude  fading  was  found  to  be  due  to  the  moon's  libration  and  was  a 

function  of  the  position  of  the  moon  in  it:  orbit.  By  comparing  the  slow  signal 

amplitude  fading  on  two  close  spaced  frequencies,  120  0  megacycles  ±0.6  per 

cent,  which  was  due  to  the  rotation  of  the  plane  of  polarization  by  the 
13  lk 

ionosphere,  Evans  9  was  able  to  deduce  the  ionospheric  electron  content  during 
periods  when  the  moon  crossed  the  meridian  at  Manchester,  England, 

Bauer  and  Daniels^  performed  similar  measurements  at  a  frequency  of 
151*11  megacycles  employing  a  bistatic  radar  system  operating  between  Belmar, 

New  Jersey  and  Urbana,  Illinois  The  absolute  values  of  the  total  electron 
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content  were  evaluated  from  these  measurements 

The  lunar  observations  cf  Hill  and  Dyce^  at  Menlo  Park,  California  on 
106.1  megacycles,  however,  did  ret.  permit  an  accurate  determination  of  the 
electron  content  because  of  the  ambiguity  in  resolving  the  number  of  polarization 
rotations  of  the  wave 

B  STATEMENT  OF  PROBLEM 

The  analytical  results  of  the  radar-lunar  measurements  conducted  at  the 
U-  S  Air  Force  Trinidad  Test  Site,  WIF,  ( 10 1-7  T°  N,  6l  6°  W)  during  the  period 
between  January  and  J lly  i960  arA  herein  reported  The  objectives  of  the 
experimental  program  were  to  determine  (ly  the  diurnal  and  also  seasonal 
variation  of  the  electron  content  in  the  ionosphere  (2)  the  reflectivity 
characteristics  of  the  me  Vs  surface 

Transmissions  were  made  at  a  frequency  of  approximately  425  megacycles 
utilizing  linear  (horizontal)  polarization  while  the  reflected  lunar  signals 
were  simultaneously  received  on  both  tne  transmitted  and  ortnogonal  polarizations. 

Th^  moon  was  automatically  tracked  from  moonrise  to  moon set  on  two  of 
the  experimental  measurements  while  partial  lunar  orbits  were  observed  on  four 
other  occasions  The  date  of  observation  of  one  of  the  partial  lunar  orbits 
corresponded  to  the  time  when  the  moon  was  at  its  maximum  southern  declination. 

The  remaining  five  observations  were  conducted  during  the  tire  when  the  moon 
was  its  maximum  northern  declination 


3 


II.  THE  DETERMINATION  OF  THE  ELECTRON  CONTENT  IN  THE  IONOSPHERE 


A  INTRODUCTION 


Prior  to  the  advent  of  rockets ,  space  vehicles  and  high-powered  radars, 
the  de termination  of  the  characteristics  of  the  ionosphere  such  as  the  distribution 
of  electron  density  with  height  and  the  electron  content  was  accomplished  by  vertical 

17 

incidence  pulse  sounders  and  radio  star  methods.  1  With  the  availability  of  these 
new  techniques,  the  state  of  knowledge  of  the  structure  and  composition  of  the 
ionosphere  has  been  greatly  enhanced. 


Electron  density  profiles  above  the  maximum  ionization  of  the  F -layer  are 

.18,  19,  20, 


readily  obtainable  from  radio  wave  transmissions  emanating  from  rockets 


21,  22, 


and  satellites 


23,  2h 


25  26 

Trom  r-f  probe -satellite  measurements  ;  and  by 

27  28  29 

means  of  incoherent  scattering  of  radio  waves  by  free  electrons.  ; 


Radio  signals  from  satellites  undergoing  doppler  frequency  changes, 


30,31 


and  Faraday  rotation  effects 


-3. 


32 


can  be  analyzed  to  yield  the  electron  content 


in  the  ionosphere.  In  addition,  the  ionospheric  electron  eontent  can  also  be  deduced 
from  radar  reflections  from  the  moon  which  display  Faraday  rotation  fades . ^ 

16,  33 


In  this  section,  the  theory  of  analysis  of  radar-lunar  amplitude  data  for 
determining  the  electron  content  in  the  intervening  medium  is  briefly  discussed. 


B .  FARADAY  ROTATION 


When  a  linearly  polarized  electromagnetic  wave  is  propagated  through  an 
ionized  medium  immersed  in  a  magnetic  field,  the  plane  of  polarization  undergoes  a 
rotation.  This  phenomenon  eommonly  referred  to  as  the  Faraday  effect  is  known  to 
take  place  for  radio  waves  traversing  the  ionosphere . 

It  can  be  shown  (see  Appendix  A  that  the  amount  of  rotation  experienced 
by  a  linearly  polarized  wave  traversing  a  two-way  path  in  the  ionosphere  can  be 
represented  by  the  function 


0  = 


L.7233  x  10 
„  2 


o 


f (h  )  H  eos  0  N  dh 


a) 


u 


where  0  is  the  rotation  in  radians,  f  is  the  transmission  frequency  in  cps,  H  is 

o 

the  magnetic  field  intensity  in  gauss ,  N  is  the  electron  density  in  electrons  per 
cubic  centimeter ,  dh  is  the  height  differential  in  centimeters ,  and  0  is  the  pro¬ 
pagation  angle;  i.e.,  the  angle  between  the  direction  of  the  earth's  magnetic  lines 
of  force  and  the  direction  of  electromagnetic  propagation. 

The  function;  f(h),  which  is  basically  the  secant  of  the  angle  between 
the  ray  path  and  the  zenith;  is  given  by 


(2) 


where  r  is  the  radius  of  the  earth  and  E  is  the  elevation  angle  of  the  antenna 
o 

beam . 

Equation  (l )  is  representative  of  the  quasi -longitudinal  mode  of  propagation. 
For  frequencies  in  the  4?5-mc  range,  this  relationship  is  valid  for  0°  <  0  <  89°  • 

It  should  be  noted  that,  for  transmissions  to  the  moon;  the  total  angular 
rotation,  0,  is  ambigous  in  that 

0  *  mt  ±  A0  (3) 

vhere  n  is  a  positive  integer.  The  parameter,  £0,  is  the  acute  polarization  angle 
normally  indicated  in  a  radar-lunar  measurement  and  is  defined  by 

=  Tan'1  —  (4) 

L  AT  J 


where  A^  and  are  the  amplitudes  of  the  signals  received  on  the  orthogonal  and 
transmitted  polarization  channels,  respectively. 

One  method  for  resolving  the  ambiguity  of  the  complete  number  of  polar¬ 
ization  rotations  is  to  employ  two  closely  spaced  frequencies.  From  the  lunar- 
echo  amplitude  fading  patterns  recorded  at  the  two  frequencies  it  is  possible  to 

infer  the  electron  content  in  the  ionosphere.  This  technique  was  employed  by 
13,  14 

Evans . 
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Since  only  one  frequency  was  available  in  the  Trinidad  radar-lunar 
experiment,  it  was  necessary  to  consider  various  methods  for  resolving  the 
ambiguity  problem.  It  was  shown  in  a  preliminary  analysis  of  the  experimental 
data  that  it  was  possible-,  by  a  process  of  elimination,  to  arrive  at  a  somewhat 

33 

reasonable  solution  to  the  polarization  ambiguousness  This  was  based  on  the 
correlation  of  the  ionospheric  electron  content,  deduced  from  the  radar-lunar 
amplitude  data,  with  the  electron  content  up  to  the  F-layer  maximum  derived  from 
vertical  incidence  ionospheric  soundings. 

The  analysis  of  the  data  presented  in  this  study  has  been  modified  to 

33 

some  extent  from  the  earlier  published  results.  That  is,  instead  of 
indirectly  arriving  at  a  solution  to  the  ambiguity  problem  by  a  process  of 
elimination,  theoretical  estimates  are  made  of  the  expected  angular  rotation 
along  various  earth -moon  paths.  These  data,  in  turn,  are  compared  with  the 
experimental  observations.  The  measured  angluar  rotation  as  defined  by  Equation 
(3)  which  best  correlates  with  the  theoretical  calculation  is  then  selected  for 
the  parameter  to  be  used  in  evaluating  the  electron  content  along  the  ray  path. 

For  all  lunar  orbits  traversing  paths  north  of  Trinidad,  1  e.,  moonrise 
occurring  at  an  azimuth  angle  less  than  90  degrees  and  moonset  at  an  azimuth 
angle  greater  than  270  degrees,  the  electron  density  profiles  evaluated  from 
vertical  incidence  ionospheric  sweep  frequency  soundings  recorded  at  the 
National  Bureau  of  Standards1  Puerto  Rico  station  were  employed  in  the  computation 
of  the  theoretical  Faraday  rotation.  The  electron  densities  deduced  from  the 
National  Bureau  of  Standards  ionospheric  sounder  station,  Bogota,  Colombia,  were 
used  in  the  analysis  of  the  lunar  southerly  orbital  path,  i.e  ,  moonrise  occuring 
at  an  azimuth  angle  greater  than  90  degrees  and  moonset  at  an  azimuth  angle  less 
than  270  degrees.  The  loci  of  points  of  the  lunar  orbits,  as  projected  on  the 
earth’s  surface  from  100-km  and  500 -km  altitudes,  are  illustrated  in  Figure  1. 

The  virtual  height  versus  frequency  records  of  the  Puerto  Rico  ionospheric 

sounder  were  converted  into  true  height  versus  electron  density  profiles  by 

3U 

means  of  the  Budden  matrix  technique  ,  the  coefficients  of  the  matrix  being 
furnished  by  E  R.  Schmerling  of  the  Ionosphere  Research  Laborabory,  the 
Pennsylvania  State  University.  The  ionosonde  data  recorded  at  Bogota  were 
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Figure  1.  Locus  of  Points  of  the  Moon's  Oroit  at  Maximum  Declination  as  Viewed  from  Trinidad 


reduced  to  true  height-electron  density  profiles  by  J.W.  Wright  cf  the  Central 
Radio  Propagation  Laboratory,  National  Bureau  of  Standards. 


Since  sweep  frequency  sounders  are  capable  of  only  yielding  the 
characteristics  of  the  ionosphere  up  to  the  height  of  maximum  ionization  of  the 
F-layer,  it  was  necessary  to  postulate  a  model  above  this  region  for  the  Faraday 
calculations.  As  a  first  approximation,  it  was  assumed  that,  above  the  F-layer 
maximum,  the  distribution  of  electron  density  with  height  (under  equilibrium  conditions 
followed  a  Chapman  model  expressed  by  the  function 


N  =  N  e 
m 


1/2 


(5) 


where  N  is  the  electron  density  at  the  level  of  maximum  ionization  and  Z  is  the 
m 

normalized  height. 

The  normalized  height  is  given  by 
h  -  h 

z  =  - £ -  (6) 

H' 

where  h^  is  the  height  of  maximum  ionization.  The  scale  height,  H',  which  is  the 
height  of  the  homogeneous  atmosphere  at  a  given  temperature,  is  defined  by 

H »  =  KT  (7) 

where  K  is  Boltzmann's  constant  (I.38  x  10  erg/degree  Kelvin);  T,  the 

-27 

temperature  in  degrees  Kelvin;  m,  the  mean  molecular  mass  of  air  (4.8  x  10  J  gm); 

2 

and  g,  the  gravitational  constant  (980.2  cm/sec  ). 

It  should  be  noted  that  the  scale  height  associated  with  the  Chapman 
expression,  Equation  (5),  is  that  of  the  neutral  particles.  In  this  study,  a 
constant  scale  height  of  80  kilometers  was  assumed  for  all  the  electron  density 
distributions  above  the  F  region. 

According  to  Equation  (1),  the  magnitude  of  the  angular  rotation 
encountered  by  a  linear  polarized  wave  transmitted  through  the  ionsophere  is 
proportional  to  H  Cos  Q  where  H  is  the  magnetic  field  intensity  and  0  is  the 
propagation  angle . 
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It  can  be 
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shown  that 


0  =  Cos 


Cos  a  Sin  1  Sin  e  Cos  I  Cos  (  y  -  D) 


(3) 


where  I  and  D  are  the  magnetic  inclination  and  declination  angles,  respectively. 
These  parameters  specify  the  direction  of  the  total  magnetic  intensity  vector  in 
space . 


The  angle,  7,  is  the  geographic  azimuth  bearing  of  the  observer's 
location  measured  with  respect  to  the  subionospheric  point,  i.e.,  the  location 
on  the  earth’s  surface  directly  underneath  the  particular  point  in  space.  The 
angle,  e,  is  simply  given  by 


e 


Sin 


-1 


r 

o 


r  +  h 

o 


Cos  E 

J 


(9) 


where  h  is  the  height  of  the  point  in  space  above  the  earth's  surface;  rQ  and  E 
are  previously  defined  in  Equation  (2). 

The  calculations  of  H  and  0,  as  reported  in  the  preliminary  results  of 
the  Trinidad  data*"'",  were  based  on  the  assumption  that  the  earth’s  magnetic  field 
may  be  approximated  by  a  centered  magnetic  dipole  having  a  geomagnetic  pole  at 
78.6°  H  and  70.1°  W. 


The  technique  employed  in  evaluating  H  and  0  in  this  report  is  as 
follows:  the  parameters,  H,  I  and  D  have  been  scaled  in  2.5-degree  and  5 -0-degree 
latitude  ana  longitude  steps  over  the  whole  earth’s  surface  from  isomagnetic  maps 
issued  by  the  U . S .  Navy  Hydrographic  Office  and  the  Canadian  Department  of  Mines 
ana  Technical  Surveys.  The  ground-observed  magnetic  data,  scaled  from  these  maps, 
are  stored  in  matrix  form  In  an  IBM-7090  digital  computer.  Linear  interpolation 
is  used  for  obtaining  the  values  of  the  magnetic  field  elements  at  other 
geographic  locations . 

When  traversing  a  path  in  the  ionosphere,  the  angles  I  and  D  at  the 
subionospheric  point  are  assumed  to  be  invariant  at  all  heights.  The  magnetic 
field  intensity  at  the  point  in  space,  on  the  other  hand,  is  derived  from  the 
surface  value  which  is  assumed  to  decrease  inversely  as  the  cube  of  the  distance 
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1 


from  the  earth's  center 


i 


it  is  believed  that  this  method  for  determining  H  and  0  is  far  more 
accurate  than  that  obtained  by  asstiming  a  dipole  model  for  the  earth’s  magnetic 
field  The  evaluation  of  the  magnetic  field  elements  based  on  spherical 

3  6 , 3? 

harmonic  expansion  is  another  method  which  could  be  employed.  However, 

utilization  of  such  a  sophisticated  technique  is  believed  not  warranted  in  the 
Faraday  rotation  problem. 


C.  ELECTRON  CONTENT  DETERMINATION  BY  RADAR -LUNAR  TECHNIQUE 


Once  the  ambiguity  of  the  total  number  of  rotations  defined  by  Equation 
(3)  is  resolved,  it  is  then  readily  possible  to  determine  the  total  electron 
content  in  the  ionosphere. 

For  transmissions  to  the  moon,  Equation  (1)  can  be  modified  to  the  form. 


0  = 


4.7233  x 


m 


/  f ( h )  H  Cos  0  N  dh 
d 


(10) 


'  CO  00 


f(h)  H  Cos  G 

r 

1  M  dh 

- 

h  Jh  J 

in  m 


where  h  ^3  the  height  of  the  F -layer  maximum  and 
m - - 

f(h)  H  Cos  0 

m 


9 

is  the  average  value  of  the  function,  f(h)H  Cos  0,  between  the  limits  of  the 
neight  of  the  F -layer  maximum  and  the  distance  to  the  moon.  As  a  first 
approximation,  it  is  valid  to  assume  that 


f(h)  H  Cos 


f (h  )  H  Cos  e 


1000  km 

h 

m 


10 


since  the  contribution  of  the  integrated  electron  density  to  the  angular  rotation 
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of  the  wave  is  considered  to  be  negligible  above  a  height  of  1GOG  km. 

The  value  of  f(h)  II  Cos  Q  as  viewed  from  Trinidad  is  plotted  as  a 
function  of  height  in  Figure  2  for  various  azimuth-elevation  orientations  of  the 
northern  and  southern  lunar  orbits.  It  is  evident  that,  at  low  angles  of 
elevation,  the  function  can  vary  over  quite  a  wide  range.  As  the  elevation  angle 
is  decreased,  there  is  a  noticeable  reduction  in  the  spread  of  its  value  as  a 
function  of  height. 

Because  of  the  extreme  variation  in  magnitude  at  low  elevation  angles, 
it  is  necessary  to  retain  f(h)li  Cos  0  within  the  integral  up  to  the  F-layer 
maximum  height.  Above  this  height,  the  function  can  be  approximated  by  an 
average  value  taken  between  the  F-layer  maximum  and  the  1000  km  height  limits. 

The  electron  content  above  the  F-layer  maximum  can  be  expressed  in  terms 
of  the  electron  content  below  the  F-layer  maximum  by  the  constant  K  where 


h  o 


rr 


(12) 
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According  to  Hisbet  "  and  Nisbet  and  Bowhill  ,  a  description  of 
the  ionosphere  in  terms  of  the  constant  K  is  not  ideal  since  the  electron 
content  below  the  F-layer  maximum,  determined  from  ionograms  may  be  in  error  due 
to  an  underestimation  of  20  to  ho  km  of  the  crue  height  of  the  F-layer  maximum. 
The  existence  of  this  error  which  is  prevalent  during  the  night  was  disclosed 
from  rocket  measurements .  The  daytime  observations  of  electron  density  profiles 
by  Ilisbet  were  in  good  agreement  with  those  obtained  from  the  true  -  height 
scaling  of  ionograms . 

A0 

In  concurrence  with  Evars  and  Taylor  ,  the  use  of  the  ratio  K  will  be 
maintained  in  this  study  because  of  the  lack  of  a  substitute  parameter  to 
describe  the  variations  of  electron  content  in  the  ionosphere. 

Foss“^  has  proposed  that  the  simplest  parameter  that  can  be  readily 
determined  for  describing  electron  content  measurements  is  the  equivalent  slab 
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Figure  2.  The  Function,  f(h)  II  Cos  9,  as  Viewed  from  Trinidad 


thickness  of  the*  total  lonosi)here,  t,  given  by 


t  N  -  / 
m  J 


h  , 

p  .ii 


H  rill  +  J  H  dh 


in 


(13) 


where  ;;  is  the  maximum  electron  density  of  the  7-layer .  The  first  integral  is 
the  electron  content  below  the  peal:  of  the  7-layer  evaluated  from  vertical 
incident*  sounders  while  the  second  is  the  electron  content  deduced  from 
Faraday  rotation  me  a*.  lements  (Equation  (10)).  It  should  be  noted  that  r  is 
essentially  the  equivalent  slab  thickness  of  the  F-layer. 
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Wright  has  shown  that, 
1  no,  W_r,  ,i  of  the  'hnrrnr. 
be  ox  ?  1  by 


if  the  electron  density  distribution  in  the 
tr  total  integrated  electron  content  can  then 


00 


This  relationship  can  be  used  to  calculate  the  scale  height,  H’,  from  tl  5 
experimental  data. 

It  is  seen  from  Equations  (13)  and  (lb)  that,  for  an  assumed  Chapman 
distribution  of  electrons,  the  equivalent  slab  thickness  is  simply  related  to  the 
scale  height  by  a  constant  factor  cf  b.133. 
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Ill  THE  CHARACTERISTICS  OF  RADIO  WAVES  REFLECTED  FROM  'THE  MOON 


A  INTRODUCTION 


The  characteristics  of  the  noon  ^s  a  radar  target  are  determined  almost 
entirely  by  the  surface  structure  of  the  moon  Although  large  scale 
irregularities  rich  as  mountains  and  valleys  are  known  to  ^xist  from  optical 
observations  the  mall  scale  irregularl ties  having  dimensions  of  the  order  of 
j.  wuve  length  ;  u.tly  mfluen  €  t re  reflection  oi  radio  waves  from  its  surface. 


In  S-dditton  irfare  r<  archness,  the  relative  velocity  of  the  moon  with 
respect  the  eartn  ml  lunar  libration  arc  also  important  factors  which 
contribute  ro  tin  characteristics  of  the  resultant  echo-  The  results  of  such 
effects  are  evident  by  variations  m  pulse  amplitude  fading,  pulse  broadening,  and 
doppler  frequency  shift. 


In  this  section,  reflection  laws  governing 
from  various  surfaces  as  determined  by  optical  means 
reflected  by  the  moon  The  probability  distribution 
lunar  echoes  urd  the  dopj^r  frequency  spread  due  to 
briefly  liseussed 

B  LUNAR  REFLECTION  LAWS 

A::ord;ng  to  the  radar  equation,  the  power 
reflect#  radiation  isotropical  1 y,  Is  given  by,  for  a 

p  -v  2  r 

i  t  w  K  C 

K  - - - 

r 

( hr  r  R  * 


the  scattering  of  radiation 
are  applied  to  radio  signals 
of  the  amplitude  of  radar - 
the  moon's  libration  are  also 


returned  from  a  target,  which 
monostatic  system, 


(15) 


where  is  "he  reccL*  i  p  the  transmitted  peak,  power;  G,  the  gain  of  the 

antenna  relative  *c  an  i 5  tro pic  antenna;  P.  the  radar-target  range;  X,  the 
transmission  wavelength;  and,  r ,  the  target  cross  sectional  area 


For  an  idealized  sphenca.  target  whose  dimensions  are  large  compared 
with  wavelengtn  and  whose  surface  is  smooth  and  perfectly  reflecting,  its 
cross  section  is  merely  "he  geometric  projected  area  given  by 

n  «•  (16) 


lU 


where  r  is  the  radius  of  the  sphere. 

When  the  target  is  not  a  totally  reflecting  body,  the  cross  section  is 
modified  to 

cs  -  It  r‘  p  (17) 

where  p  is  the  power  reflection  coefficient  of  the  surface. 

It  is  evident  from  the  Fresnel  diffraction  theory  that  the  amplitude  of 
a  spherical  wave  returned  from  a  smooth  body  is  equal  to  one-half  that  due  to 
the  first  Fresnel  zone  acting  alone.  This  applies  to  an  aperture  with  an  infinite 
number  of  zones.  Since  the  first  Fresnel  zone  has  a  radius  of  (  R  X/2)  '  ,  in 

the  case  of  the  moon  (if  it  were  a  smooth  body),  only  a  small  portion  of  the 
moon’s  surface  would  contribute  to  the  received  reflected  energy.  The  first 
Fresnel  zone  occupies  a  portion  of  the  moon  nearest  to  the  observer  so  that  a 
pulse  reflected  from  the  moon  would  not  be  extended  by  reflections  from  the  rear 
portions  of  the  moon,  in  general,  it  can  be  stated  that  a  pulsed  electromagnetic 
wave  incident  on  a  smooth  sphere  is  reflected  without  any  distortion,  i.e.,  pulse 
width  of  the  reflecte  1  wave  is  identical  to  that  of  the  incident  wave. 

With  regard  to  a  spherical  target  whose  surface  is  comparatively  rough, 
the  equivalent  echoing  area,  o^,  is  given  by 

2 

°r  -  *  r  p  q  (18) 

where  q  is  the  directivity  of  the  target.  The  magnitude  of  the  directivity 
factor  is  dependent  upon  the  scattering  law  applicable  to  the  rough  surface. 

Although  the  radiation  is  scattered  in  all  directions,  there  is  a  greater 
increase  in  the  amount  of  energy  scattered  back  in  the  echoing  direction  for  a 
rough  sphere  than  for  a  smooth  sphere  of  .he  same  geometric  size. 

The  two  scattering  laws  considered  in  this  report  are  the  Lambert  and 
the  Lommel-Seeliger  lav. 


15 


The  Lambert  law,  which  applies  to  the  scattering  from  a  diffuse 
surfaee  having  irregularities  on  the  order  of  a  wavelength,  states  that  the 
scattered  energy  in  any  direction  is  proportional  to  the  cosine  of  the  angle 
between  the  incident  ray  and  the  normal  to  the  surface  and  to  the  cosine  of  the 
angle  between  the  scattered  ray  and  the  normal. 

Tn  applying  the  Lambert  law  to  the  scattering  from  a  rough  sphere,  it 
is  shown  in  Appendix  B  that,  for  an  incident  pulse  of  finite  length,  the  received 
reflected  power  at  any  instant  of  time  is  merely  the  sum  of  all  the  energy 
contributions  from  a  finite  band  around  the  sphere.  The  received  power,  assuming 
reflection  back  in  the  direction  of  the  radiating  source,  can  be  represented  by 
the  function 


P  = 
r 


9  9 

8  P  C  x  a 
_ t _ s 

(Uk)3  R4 


l  31 


Sin  0  Cos"  9  d0 


(19) 


9. 


where  Q  is  the  angle  measured  from  the  radial  line  passing  through  the  point  of 
tangeney  containing  the  incident  wave  front  and  the  sphere.  The  limits  of 
integration  are,  for  the  condition  t  >  T 

4"  (  t-T) 


r  - 


0 i  =  Cos 


-1 


(20) 


ana 


0  =  Cos 


-1 


r  r  - 


ct 


(21) 


where  T  is  the  pulse  width  and  t  is  the  travel  time  of  the  leading  edge  of  the 
pulse  along  the  moon’s  surface  measured  from  the  point  of  tangeney.  When 
T  >  t ,  -  C  while  0^  remains  unchanged. 

It  is  evident  from  Equations  (15)  and  (19)  that,  for  a  pulse  width 
greater  than  11.6  milliseconds,  i.e.,  time  required  to  travel  a  distance  of  twice 
the  radius  of  the  moon,  the  signal  power  received  from  a  rough-surfaced  moon, 
obeying  Lambert's  law,  would  be  8/3  times  greater  than  that  from  a  smooth  moon. 
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It  t.h  r  \k  >_a  ft  1  'ivity  fuc*or,  ic scribed  in  Equation  (10 ),  for  the 
Lumber:  -are,  t  equal  igo/M 

It  Lomme i -See l.Lger  law  refers  to  scattering  taking  place  from  a  rough 
urf  ■  navmg  irregularities  tJau~  arc  large  compared  with  wavelength.  The 
cue- r g  /  i  ;  it  t  ~re  1  from  al  1  regi'  nft  of  t  ri riiri  aee  the  same.  Thus,  when 
'  cr. ide  r  In;;  ruuio  wave*  re  flee  tl  n  l  rom  su:*h  a  surface,  the  received  power  is  only 
proport lonai  tv  tin  :c  me  cr  rre  angle  between  the  incident  ray  and  the  normal 
r i  r  c  z  ir  fae  c  . 

Av  ..h  w:  m  Appendix  L  ♦hr  r:  rved  power  can  then  be  expressed  by 


- 

(  *• n  i  R 

I-  f  ilews .  there:'  rc  Thn1  for  ouiSe  lengths  greater  than  11.6 
auMlser  :nd:  .  tvim  the  power  is  received  from  a  moon,  approximating  Lommel- 
3  c  i  :g*.  i  ea:  •  ring  ~h&r.  from  t  mot  ih  moon. 

Wtf  .  tgard  to  the  effect  of  scattering  from  a  rough  surface  on  the 
>  .  “e  miv  ill  portions  of  a  rough  body  contribute  to  the  total 

ref.1  ter?  signal  the  returned  euc  woteJUl  be  extended  in  length,  i.e.,  the  pulse 
lengthening  n  it:  ♦.  ]nal  tc  the  time  it  takes  the  wave  to  travel  from  the  nearest 
Surface  t;  *he  limb  and  then  beck  again  to  the  surface.  Since  the  radio-depth 
i  the  moon  !  -v  M  .1,  mil  1  i  ecords,  a  pulse  o:  length,  T,  incident  on  a  rough 
m -r  .  wc  no  be  ?  ongated  to  a  length,  T  -  11. h  milliseconds,  after  reflection. 

STATIST! tAL  DISTRIBUTION  OF  THE  AMPLITUDE  OF  LUIIAF  ECHOES 

1  •  Farid  m  Fading 

T;  ^  rapid  fluctuations  of  rauar  .  uises  reflected  from  the  moon's 
irfaie  nr?  usually  attributed  to  the  li  brat  ion  which  is  defined  as  the 
-  pil:ut  ry  nr.ion  cl  ~v  e  r  oori  a  bom  m  axic  which  itself  changes  with  time.  There 
r* re  M  prm1:  ,  causes  of  libration  (‘1  «  libration  in  latitude:  this 
r*.  j'<  from  the-  ih  *  mat  the  sxu  if  rotation  of  the  moon  is  tilted  about 


f  r' 

/  2  in  6  Ccs  0  dv  (22) 

A 

1 


1  ^ 


6.5  degree?  to  the  plane  uf  revolution  Thus,  from  a  fixed  position  on  the 
earth  varying  amounts  of  the  moon’s  northern  Hemisphere  and  southern  hemisphere 
are  \isible;  (2)  iibration  in  longitude;  this  is  due  to  the  fact  that  the 
period  of  rotation  of  the  moon  is  constant,  but  the  speed  along  the  orbit  varies 
In  accordance  with  Kepler's  second  law  fcr  planets;  (3)  diurnal  Iibration;  this 
is  the  result-  of  the  rotation  ot  the  earth.  As  the  earth  rotates,  the  position 
f  the  moo a  relative  to  an  observer  on  the  earth  is  continually  changing. 

•The  effe  r  f  the  moon’s  Iibration  on  the  reflected  signal  can  be 
u  In  mod  by  is  sum  mg  that  the  moon’s  surface  consists  of  a  random  number  of 
scat terers .  V he  amplitude  of  the  reflected  pulse,  as  observed  on  the  earth’s 
^urface,  m  the  resultant  of  the  signals  reflected  from  each  of  the  scattering 
element^.  3m  e  the  moon  undergoes  ar:  apparent  rocking  motion, or  Iibration, 
the  signal?  scattered  from  various  jarts  of  its  surface  are  continuously 
undergo  mg  ransom  changes  in  phase  and  amplitude  which,  in  turn,  produce 
f ;uc "u&t  t ens  m  the  resultant  signal. 


For  ur; ace  irregularities  which  scatter  incident  radiation  with  random 
phases  ml  sjn:ii*udes,  rne  probability  ot  occurrence  any  resultant  signal  is 
therefore  giver  by  me  Rayleigh  probability  distribution  law.  The  concept  of  the 
Rayleigh  distribution  m  maintained  provided  that  the  number  of  scattering 
area:  1  large  .  .e.,  at  least  on  the  order  of  10  or  greater. 


The  probability,  FfR  dR ,  of  finding  an  amplitude  between  R  and  (R  t-  cLR) 
fnerefore  given  by 


Pi  R 


~  4 


(23) 


wb:  'e  R  is  the  amplitude  of  the  resultant  scattered  signal  at  any  instant  of  tire 
and  $  n  a  constant  iefmed  m  terms  of  the  first  or  second  moments  of  R  by 


or 


n  ^  12 

-7T1—) 


2 


-  2* 


(2U) 

(25) 
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It  can  be  readily  shown  that  t  can  also  be  expressed  by 


Rm  -  (t)l/2  (26) 

where  R  is  the  most  probable  value  of  R,  i.e.,  the  value  of  R  when  P(R)  is  a 
maximum . 


It  should  be  noted  that  the  distribution  of  the  amplitude  of  random 

2*3 

noise  emerging  from  a  bandpass  filter,  as  discussed  by  Rice  ,  follows  the 
Rayleigh  law. 

2 .  Random  Fading  Contaminated  with  a  Steady  Signal 

The  condition  in  which  a  steady  signal  is  superimposed  on  a  signal 
which  is  the  resultant  of  elementary  contributions  originating  from  random 
scattering  areas  is  of  considerable  interest  since  it  approximates  a  possible 
situation  prevailing  in  the  lunar  reflection  of  radio  waves. 

The  steady  signal  corresponds  to  a  specularly  reflected  wave  coming 
from  the  first  Fresnel  zone  or  from  a  relatively  small  number  of  smooth  surface 
areas  in  the  region  of  the  first  Fresnel  zone.  The  random  signal  is  caused  by  the 
lunar  libration  or  by  random  surface  irregularities.  If  the  number  of  smooth 
surface  areas  becomes  very  large,  it  is  possible  for  the  signals  reflected  by 
these  areas  to  undergo  cancellation  and  reinforcement.  Thus,  instead  of  specular 
reflection  taking  place,  a  random  type  noise  could  result. 


The  probability  density  distribution  function  that  describes  the 
envelope  of  resultant  signal  can  be  written  as,  according  to  Rice's  theory 
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of  random  noise  , 


r  V  v  ) 


V  B 


(27) 


where  B  is  the  amplitude  of  the  steady  signal;  V,  the  envelope  of  the  resultant 
amplitude  comprised  of  the  steady  signal  (B)  and  the  random  signal  (R); 
constant  defined  by  Equations  (2h),  (25),  or  (26). 

The  parameter,  I  (VB/V),  is  the  modified  Bessel  function  of  the  first 
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kind  of  zero  order,  given  by 


i  (-?■)- 

o'  \Jr  ' 


,  VB  \  2n 
^~2 —> 


(28) 


n=0 


During  B  =  b-y/\jT,  it  follows  than,  when  b=0,  Equation  (27) 
simplifies  to  the  Rayleigh  distribution  law. 


When  b>3,  the  Bessel  function  can  be  replaced  by  its  asymptotic 


expression 


I.  (-gr-)  -  (- 


"o  'VV  '  '  2n  b  V 

and  then  Equation  (27)  simplifies  to 


a/2  _ b  V 

)  e  ^ 


(29) 


-\A' 


P(V)*=  -i-  e--(V-bV^r 

VSfT  ^  b.,y  e 


(30) 


_  l/2 

It  is  evident  that,  apart  from  the  term  (v/b  V^)  this  distribution  is 
Gaussian  with  a  standard  deviation  of-/t*  When  b  is  equal  to  or  larger  than  two, 
the  presence  of  a  steady  signal  is  indicated. 

When  applying  Equation  (27)  to  experimental  data,  it  is  necessary  to 
obtain  the  value  of  Jr.  If  the  data  to  be  analyzed  is  in  amplitude  form, 
containing  both  steady  and  random  signals,  then  according  to  Equation  (24), 

\Jr  may  be  approximated  by 


♦ 


v  2 


(31) 


It  follows  from  Fquation  (25)  that  if  the  data  is  in  the  form  of 
cross  section,  a,  or  power, 


(32) 


It  should  be  mentioned  that  the  validity  of  the  approximations  for 
o  taining  ;  from  experimental  data,  Equations  (31)  and  (32),  is  improved  for 
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b  <  2. 


In  applying  the  theory  of  random  noise  to  the  study  of  the 
characteristics  of  the  lunar  surface,  it  would  appear  that,  if  the  moon  were 
completely  smooth,  b  would  be  on  the  order  of  three  or  greater  and  if  it  were  a 
rough  body,  b  would  most  likely  be  zero.  Assuming  that  the  moon  contained  both 
rough  and  smooth  regions,  it  could  be  expected  that  b  would  lie  between  zero  and 
three.  However,  it  should  be  mentioned  that,  for  this  particular  situation,  it 
would  be  difficult  to  determine  the  extent  of  the  actual  smooth  and  rough  areas. 

Specular  reflection  can  contribute  to  the  resultant  signal  when  the 
surface  normal  is  in  the  direction  of  the  observation  site;  such  areas  are  most 
likely  contained  on  a  cap  of  the  moon  closest  to  the  earth.  The  specular 
returns  can  be  highly  perturbed  either  by  rough  scattering  areas  located  in  the 
near  vicinity  of  smooth  regions  or  by  the  presence  of  a  large  number  of  isolated 
smooth  surfaces. 

Another  factor  which  must  be  considered  in  the  analysis  of  lunar 
amplitude  data  is  the  time  interval  over  which  the  probability  distribution  is 
applied.  Since  the  lunar  libration  affects  the  reflected  signal  amplitude,  the 
time  over  which  any  statistical  analysis  is  applied  should  be  small  compared  to 
•the  time  variation  of  the  libration  rate. 

In  short,  an  evaluation  of  the  probability  density  distribution  of 
lunar  echo  amplitudes  may  not  result  in  a  true  indication  of  the  surface 
characteristics  of  the  moon.  Other  factors  such  as  pulse  lengthening  and  the 
characteristics  of  the  amplitude  decay  of  the  lengthened  pulse  should  be  examined 
concurrently. 

D.  DOPPLER  EFFECTS  OF  LUNAR  ECHOES 

Radio  waves  reflected  from  the  moon  are  found  to  exhibit  both  a 
doppler  shift  and  a  doppler  spread  in  the  transmission  frequency. 

The  doppler  shift  is  brought  about  by  the  fact  that  the  orbit  of  the  moon 
around  the  earth  is  elliptical.  Thus,  there  is  a  component  of  the  moon’s  orbital 
velocity  directed  toward  the  center  of  the  earth.  Since  the  earth  is  also 
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rotating,  there  is  an  additional  velocity  component  along  the  earth -moon 
line-of -sight,  its  magnitude  being  a  function  of  the  latitude  of  the  observation 
site  and  the  moon's  position  on  the  celestial  sphere.  The  resultant  of  these  two 
velocity  components  produces  the  frequency  shift  in  the  reflected  echo. 

For  a  monostatic  system,  the  doppler  shift  imposed  on  a  radar  signal 
reflected  from  the  moon's  surface  is  given  by 


2  R  f 


(33) 


where  f  is  the  difference  between  the  apparent  reflected  frequency  and  the 
transmitted  frequency,  f;  R,  the  time  rate  of  change  of  the  distance  between  the 
observation  site  and  the  surface  of  the  moon;  c,  the  free  space  velcoity. 

1+4 

According  to  Fricker  ,  the  relative  velocity,  R,  can  be  expressed  by 

Do  <Do  -  ro  C°=  ')  -Do  r„5T  (Cos  v> 

- (3M 


R  =  - — 

where,  as  shown  in  Figure 

D  = 


D 


D  +  r  -  p  r  D  Cos  j 
o  o  ^  o  o 


1/2 


(35) 


and  is  the  distance  from  the  center  of  the  earth  to  the  center  of  the  moon;  D, 

the  distance  from  the  observation  site  to  the  moon’s  center;  r  ,  the  radius  of 

o 

the  earth;  and  y,  the  angle  between  the  line  connecting  the  earth-moon  centers 
and  the  radial  line  from  the  earth  center  to  the  observation  site. 

The  terms  in  Equations (3^) and (35  )can  be  computed  from  the  following 
relationships ; 


Do  Sin  n 


D  =  -  D  Cot  n 
O  o 


d  n 


dt 


(36) 

(37) 
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EARTH 


Figure  3*  Earth-Moon  Geometry 
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d6 

dt 


(38) 


-it ( Cos  +}  = 

-  Cos  0  Cos  6  Sin  (LHA)  (LHA) 

CL  0 

-jt  (LIIA)  -  w--af  <“> 


Cos  6  Sin  0-Sin  6  Cos  0  Cos  (LHA) 


(39) 


where  II  is  the  horizontal  parallax  angle  of  the  moon,  b,  the  declination  of  the 
moon;  LHA,  the  local  hour  angle  of  the  moon,  RA,  the  right  ascension  of  the 
moon;  0,  the  latitude  of  observation  site;  and  cj,  the  angular  velocity  of  the 
earth.  The  moon's  celestial  coordinates  and  the  horizontal  parallax  are 
tabulated  in  the  Nautical  Almanac 

The  libration  of  the  moon  causes  a  variation  of  range -rate  in  the  radial 
distances  from  the  observation  site  to  points  on  the  moon’s  surface.  The 
differences  in  the  velocity  components  from  the  various  parts  of  the  moon's 
surface  give  rise  to  a  doppler  spread  in  the  reflected  signal  frequency.  In 
addition,  doppler  spreading  also  results  due  to  the  movement  of  the  observation 
site  toward  the  center  of  the  lunar  disk.  The  total  spread  is  the  resultant 
of  both  these  effects. 
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The  total  doppler  spread,  A  f^,  as  derived  oy  Fricker  ,  can  be 
given  in  terms  of  the  function 


a  ^  2  f  K  r  p 

A  f  =  -  \  JC—  Sin  v 

d  c  L  TB 


Cos  v 
TL 


f  K  r  R 
c  R1- 


(40) 


where  r  is  the  radius  of  the  moon 


the  total  effective  libration  rate 


in  latitude ; A- >  the  total  effective  libration  rate  in  longitude;  v,  the  angle 
XL 

between  a  location  on  the  perimeter  of  the  moon  and  the  moon's  equator;  and  K, 
the  fractional  radius  of  the  moon,  i.e.,  K  =  a/r  as  illustrated  in  Figure  3* 


f 


24 


The  first  term  in  this  expression  applies  to  the  doppler  spread  due 

to  the  moon's  libration  while  the  second  is  that  due  to  the  motion  of  the 

observer.  In  general,  the  latter  term  can  be  neglected  for  libration  rates 
-7 

greater  than  10  radians  per  second. 

When  the  motion  of  the  observer  is  not  taken  into  account,  the  maximum 
doppler  spread  occurs  at  an  angle,  v  ,  where 


v 

o 


Tan 


-1 


~Jy 


TB_ 

TL  -I 


(41) 


The  expressions  defining  the  parameters,  v,JL  and  a  ,  are  not 

ih  iii 

presented  in  this  report  because  of  their  complexity.  Computational  procedures 
for  evaluating  the  terras  are  contained  in  Reference  44. 


With  regard  to  the  effect  of  libration  on  lunar  echoes,  if  the  moon 
were  a  perfectly  smooth  body,  the  libration  would  produce  a  slight  doppler 
spread  within  a  Fresnel  zone,  its  magnitude  being  directly  proportional  to  the 
fractional  radius,  K,  <r  dimensions  of  the  zone  and  varying  as  a  function  of  the 
angle,  v.  Since  the  first  few  Fresnel  zones,  where  K  is  very  small,  contribute 
to  most  of  the  energy  of  the  reflected  signal,  the  doppler  spread  would  be 
expected  -o  be  insignificant.  The  amplitude  of  the  echo  would  have  a  slight 
modulation  with  a  bandwidth  determined  by  the  doppler  spread  over  the  contributing 
Fresnel  zones. 


For  a  moon  consisting  of  smooth  regions  (such  as  the  first  few  Fresnel 
zones)  and  rough  scattering  areas  (such  as  mountains),  it  may  be  possible,  by 
means  of  the  libration  effects,  to  reveal  the  existence  of  the  scatterers.  Since 
the  lunar  libration  introduces  a  doppler  spread  and  amplitude  fluctuation  in  the 
ref lected  signal,  the  magnitude  of  the  frequency  components  in  the  doppler  spread 
or  the  dominant  frequencies  in  the  power  density  spectrum  can  be  used  to 
determine  the  fractional  radius  of  the  moon  which  best  fits  the  experimental  data. 
This  parameter,  in  turn,  can  then  be  compared  with  photographs  of  the  lunar  face 
for  correlation  with  known  mountain  ranges. 

The  feasibility  of  mapping  the  surface  of  the  moon  by  means  of  doppler 
shift-power  spectrum  measurements  has  been  demonstrated  by  Pettengill. 
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IV.  EXPERIMENTAL  MEASUREMENTS 


The  experimental  observations  discussed  in  this  report  were  conducted  at 
Trinidad,  WIF  between  January  and  July  I960  with  a  high-powered  pulsed  radar 
operating  at  a  frequency  of  425  megacycles,  Linear  (horizontal)  polarization  was 
employed  on  transmission  while  the  reflected  radar  signals  were  simultaneously 
received  on  both  the  transmitted  and  the1  orthogonal  (vertical)  polarization. 

Automatic  tracking  of  the  moon  was  accomplished  by  means  of  an  orbit 
programmer  which  continuously  positioned,  in  three-minute  intervals,  an  84-foot 
diameter  steerable  antenna  so  as  to  follow  tne  theoretical  path  of  the  moon. 

The  theoretical  geometric  orientation  of  the  moon,  as  viewed  from  Trinidad, 
for  moonrise  commencing  on  12  January  L960  is  shown  in  Figure  4.  The  method  for 
determining  the  azimuth,  elevation  and  time  coordinates  of  a  celestial  body  for 
tracking  purposes  is  discussed  ir  Appendix  C 

The  radar  data  pertaining  to  the  lunar  echo  return  were  recorded  in  digital 
1'orm  on  magnetic  tape  in  order  to  facilitate  data  processing  and  analysis.  In 
addition,  photographic  recordings  w^re  made  of  the  various  combinations  of  the 
amplitude,  range  and  time  coordinates  schematically  illustrated  in  Figure  5* 

A  tracing  of  a  typical  A  scope  photograph  of  a  radar-lunar  echo  recorded  at 
Trinidad  is  shown  in  Figure  6.  The  upper  trace  is  the  horizontally  polarized 
signal  while  the  lower  is  the  vertical  It  is  noted  that  the  pulse  return  is 
stretched  out  to  approximately  12  milliseconds  in  length  and  that  the  echo  shapes 
of  both  polari  zatiors  are  comparatively  alike  The  maximum  radar  cross  section  of 
126.3  db  above  one  square  meter  is  obtained  on  horizontal  polarization  This 
value  compares  reasonably  well  with  the  geometric  projected  lunar  disk  area  of 
129-8  db  above  one  square  meter  (radius  cl  moon  =  1740  km)  Assuming  a  smooth 
moon  with  a  power  reflection  coefficient  of  0  15>  the  theoretical  echoing  area 
should  be  on  the  order  of  121  6  db  above  one  square  meter 

A  lunar  echo  showing  p>ulse  distortion  and  also  pulse  lengthening  is  contained 
in  Figure  7*  For  this  particular  example,  tie  signal  amplitude  is  only  about  ll4 
db  above  one  square  meter 
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Figure  4.  Tne  Orbit  of  the  Moon  at  Maximum  Northern  Declination  as  Viewed  from  Trinidad 


Figure  6.  A  Typical  Lunar  Echo  Recorded  at  Trinidad  at  425  Me  Shoving  Folse  Lengthen! r. 


A  selected  sample  of  an  amplitude  vs  time  film  record  of  moon  echoes 
displaying  Faraday  rotation  is  illustrate!  m  Figure  8  It  is  noted  that  the 
reflected  signal  is  predominantly  oriented  in  the  vertical  polarization  receiver 
channel  during  the  pei loi  wner  horizontal  polarization  vas  employed  on  transmis¬ 
sion. 

An  amplitude  vs  tim^  fading  record  illustrating  the  presence  of  lunar 
libraticn  is  shown  in  Figure  9  It  is  seer  that  the  envelope  of  the  back  portions 
of  the  pulse  returns  exhibits  a  much  higher  fading  rate  man  the  front  portion. 

It  is  interesting  to  note- that  the  fading  patterns  cf  the  main  echo  pulse  (the 
large  signal  amplitude )  obaer  -ed  on  both  polarizations  are  identical  The 
significance  of  this  is  that  the  polarization  of  the  incident  pulses  on  reflection 
from  the  front  part  of  the  mcc-n  is  .-laimamed*  which  implies  that  this  region  of 
the  moon  s  suriace  appears  reasonably  smooth  at  a  frequency  of  425  megacycles. 
Depolarization  of  a  signal  should  ci  :ur  on  reflection  from  a  rough  surface.  It 
would  follow  that  the  amplitude  fading  patterns  observed  on  orthogonal  polariza¬ 
tions  would  be  somewhat  different 

A  typical  range  vs  time  photograph  f  lunar  echoes,  revealing  pulse  lengthening^ 
is  shown  in  Figuie  10  The  s'lgr.a  mdi  a  ted  by  the  reduction  of  the  echo 

intensity  in  the  back  portion  cf  the  range  scale  is  mainly  due  to  the  moon’s 
libration 

The  schedule  cf  the  radar-lunar  observations  performed  at  Trinidad  is 
presented  in  Table  1  The  tota^  time  of  experimental  measurements  was  14-3/4 
hours,  of  which  less  than  1C  per  cent  (the  March  i960  data)  vere  conducted  when 
the  moon  was  at  its  maximum  southern  decimation  The  remaining  observations 
were  made  during  the  time  when  the  mccn  was  at  its  maximum  northern  declination. 
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Figure  8.  Amplitude  vs*  Time  film  Record  of  Lunar  Echoes  Recorded  at  Trinidad  at  ^25  Me 
Displaying  Faraday  Rotation,  8  February  i960,  1403  EST. 


9848 


\ 


<0  X 

O 

<0 

- 

3 

1 

n 

1-  < 

mN<Oj 

X  ISI 

-Z  e> 

g3j 

Q. 

2 

< 

UJ  — 

> (r 
< 

-1 

O 

0 

CL 

CL 

33 


Figure  9.  Amplitude  vs.  Time  Film  Record  of  Lunar  Ec.oes  Recorded  at  Trinidad 
at  1+25  Kc  Displaying  Lunar  Libration,  12  January  i960,  1617  EST. 
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Figure  10.  A  Typical  Range  vs.  Time  Photograph  of  Radar -Lunar  Echoes 
Recorded  at  Trinidad  at  U25  Me  12  Jan  lary  i960,  1  1?  3ST. 
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V.  DATA  ANALYSIS 


A.  CHARACTERISTICS  OF  THE  IONOSPHERIC  MEDIUM 
1 .  Preliminary  Analysis 

In  the  initial  planning  of  this  experiment  involving  the  moon,  in 
studying,  by  means  of  radar  techniques,  the  diurnal  variation  of  the  electron 
content  in  the  ionosphere,  it  was  reveuLed  that,  when  the  moon  was  at  its 
maximum  northern  declination,  its  orbit,  as  viewed  from  Trinidad,  intersected 
regions  in  the  ionosphere  where  the  direction  of  the  magnetic  field  lines  was 
perpendicular  to  the  direction  of  propagation  ^  Since  only  one  frequency  was 
available  for  the  experiment,  it  was  believed  that:  (l)  for  the  condition  of 
perpendicularity  with  the  magnetic  field  lines,  the  angular  rotation  of  the  plane 
of  polarization  of  a  linearly  polarized  wave  (due  to  the  Faraday  effect)  would  be 
negligible  for  transmissions  on  an  earth-moon  path;  (2)  it  would  be  quite  feasible 
to  resolve  the  ambiguity  of  the  complete  number  of  rotations  of  the  plane  of 
polarization  by  simply  counting  the  maxima  and  minima  of  the  polarization  angle 
of  the  resultant  signal. 

The  preliminary  analysis  of  the  experimental  data  indicated  that,  for 
transmissions  through  the  entire  ionosphere,  the  condition  of  completely  trans¬ 
verse  propagation  is  maintained  in  only  a  very  small  region  in  space.  The 
Faraday  rotation,  even  when  the  condition  of  perpendicularity  of  the  propagation 
direction  with  the  earth's  magnetic  field  is  satisfied,  can  still  be  appreciable. 

The  determination  of  the  electron  content  in  the  ionosphere  was, 

therefore,  accomplished  from  the  lunar-echo  polarization  data  by  an  indirect 

nn 

method.  All  possible  solutions  to  the  problem  were  considered  and,  by  a 
process  of  elimination,  based  on  the  comparison  of  the  estimated  electron  content 
with  the  electron  content  below  the  F-layer  maximum  deduced  from  ionospheric 
soundings  at  Puerto  Rico,  a  reasonable  solution  was  obtained 

The  magnetic  field  strength  and  orientation  employed  in  that  analysis 
were  calculated  on  the  assumption  that  the  earth's  magnetic  field  could  be  ap¬ 
proximated  by  a  magnetic  dipole  located  at  the  center  of  the  earth  and  having  a 
geomagnetic  pole  at  78  6°  N  and  70.1°W 
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A  re-evaluation  of  the  preliminary  work  resulted  in  a  more  direct 
method  for  solving  the  ambiguity  in  the  polarization  angle,  i  e.,  by  the 
comparison  of  the  estimated  theoretical  angular  rotation  with  the  experimental 
measurements  In  addition,  the  magnetic  dipole  model  for  the  earth's  magnetic 
field  was  replaced  by  the  surface  values 

The  lunar  orbit  of  the  moon  and  the  contours  of  perpendicularity  with 
the  earth's  magnetic  field  as  viewed  from  Trinidad  are  presented  in  Figure  11. 

It  is  seen  that  the  magnetic  field  orientations  based  on  the  dipole  model  and 
surface  data  are  slightly  different  The  latter  is  considered  to  be  the  more 
accurate  representation. 

2  Total  Electron  Content 

The  average  acute  angular  rotation  of  the  signals  reflected  from  the 
moon  during  the  observation  of  12  January  i960  is  presented  by  the  0  =  £$  and 
the  A  curve  of  Figure  12  Because  of  the  inability  of  an  antenna-receiver 
system  to  determine  the  true  orientation  of  the  incident  radiation,  it  is  neces¬ 
sary  to  examine  the  other  possible  values  of  the  polarization  angle  defined  by 
Equation  (3)  and  also  plotted  in  Figure  12 

The  intersection  of  any  two  "  0  "  curves  such  as  the  0  =  £0  and 
the  0  =  tt  -  A  0  curves  at  1800  E3T  always  occurs  when  the  polarization  of  the 
received  signal  passes  through  either  a  90- degree  or  a  0-degree  angular  rotation. 
It  should  be  noted  that  an  intersection  could  result  either  from  a  phase  reversal 
or  from  a  continuing  monotonic  phase  change  For  example,  there  are  two  possible 
paths  which  each  of  the  0  curves  could  follow;  the  0  =  a0  curve  could  have 
continued  along  the  B  path  instead  of  the  A  path 

The  theoretical  estimate  of  the  magnitude  of  the  ionospheric  Faraday 
rotation  that  could  have  occurred  at  the  various  azimuth-elevation  orientations 
of  the  lunar  orbit  is  also  shown  in  Figure  12  These  calculations  are  based  on 
the  true  height- electron  density  profiles  deduced  from  the  vertical  incidence 
ionospheric  soundings  taken  at  15-minute  intervals  at  Puerto  Rico  Above  the  F- 
layer  maximum,  the  electron  density  distribution  is  assumed  to  follow  a  Chapman 
model,  as  defined  by  Equation  (5)>  having  a  constant  scale  height  of  80  kilometers 
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MAGNETIC  DIPOLE  MODEL 


T-e  Lunar  Orbits  and  Contours  of  Perpendicularity  With 
the  Earth's  Magnetic  Field  as  Viewed  from  Trinidad 
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Figure  12.  Average  Faraday  Rotation  from  Trinidad  Radar-Lunar  Observation,  12  January  i960 


for  the  neutral  particles 


It  is  quite  evident  that  the  theoretical  calculations  correlate  to  a 
high  degree  with  the  experimenual  data  represented  by  the  0  *  n  -  £0  and  the 
A  curves.  Thus,  these  two  curves  are  chosen  as  the  most  likely  correct  solution 
to  the  ambiguity  problem.  It  is  of  interest  to  note  that  a  similar  conclusion 
was  arrived  at  in  the  early  analysis  of  the  radar-lunar  data 

The  electron  content  in  a  o.ie-square  centimeter  vertical  column 
through  the  ionosphere  is  depicted  in  Figure  13  The  experimental  plot  is  derived 
from  Equation  (10)  in  conjunction  with  the  values  of  0  specified  by  the 
0  -  jt  -A0  and  the  A  curves  It  is  seen  that,  in  the  late  afternoon,  the 

total  electron  content  determined  from  the  integration  of  the  combined  electron 
density  profile  (from  ionospheric  soundings)  and  a  Chapman  distribution  is 
slightly  higher  than  the  experimental  measurements,  while  in  the  early  evening 
the  reverse  takes  place.  This  would  imply  that,  for  this  particular  day,  the 
scale  height  during  the  daytime  should  be  less  than  80  kilometers  and  that,  during 
the  evening,  the  ionosphere  could  perhaps  be  represented  by  a  Chapman  model  with 
a  constant  scale  height  in  excess  of  80  kilometers. 

The  average  Faraday  rotation  observed  on  8  February  i960  is  shown  in 
Figure  lb-  It  is  seen  that  the  acute  polarization  angle  reaches  a  minimum  at 
1700  EST  which  would  infer  that  the  polarization  of  the  received  lunar  echo  at 
that  time  was  the  same  as  the  transmitted  polarization. 

An  examination  of  the  theoretical  angular  rotation  with  the  experi¬ 
mental  measurements  reveals  that  the  0  n  n  -  £0  and  the  A  curves  best  cor¬ 

relate  with  the  theoretical  data 

A  plot  of  the  total  electron  content  derived  from  the  radar-lunar 
obseivations  and  from  the  Puerto  Rico  ionosonde  data  is  given  in  Figure  15 •  It 
is  evident  that,  as  in  the  case  of  the  12  January  results,  a  constant  scale 
height  of  80  kilometers  above  the  F-layer  maximum,  assuming  a  Chapman  model,  is 
too  low  for  the  evening  The  daytime  data,  on  the  other  hand,  are  reasonably  well 
correlated 

The  average  angular  rotation  measured  on  19  March  i960  when  the  moon 


Figure  lb.  Average  Faraday  Rotation  from  Trinidad  Radar -Lunar  Observation,  8  February  i960 


Figure  15.  Comparison  of  the  Measured  Total  Electron  Content  in  a  One  Square  Centimeter 

Vertical  Column  Through  the  Ionophere  with  Theoretical  Estimate,  8  February  i960 


traversed  a  southerly  orbit,  with  respect  to  Trinidad,  is  shown  in  Figure  l6. 

The  theoretical  computations  are  based  on  vertical  incidence  ionospheric  soundings 
recorded  at  Bogota,  Colombia  and  an  assumed  electron  density  distribution  of  the 
Chapman  form  with  a  scale  height  of  80  kilometers  above  the  F-layer  maximum. 

The  conversion  of  the  Bogota  virtual  height  versus  frequency  records  into  true 
height  versus  electron  density  was  performed  by  the  National  Bureau  of  Standards. 

The  dashed  lines  between  05 -*5  and  06^5  hours  indicate  missing  experi¬ 
mental  data  Based  on  the  theoretical  calculations,  the  measured  angular  rota¬ 
tion  is  assumed  to  have  a  value  of  90  degrees  at  approximately  0615  hours  It  is 
quite  evident  therefore,  that  the  actual  Faraday  rotation  encountered  in  the 
ionosphere  can  be  described  by  the  0  =  A  0  ,  B  and  C1  curves 

The  variation  of  the  total  electron  content  deduced  from  an  assumed 
angular  rotation  given  by  the  0  =  A  0  ,  B  and  C  *  curves  is  contained  in  Figure 

17.  The  experimental  measurement  of  the  total  electron  content  is  found  to  be 
higher  than  that  theoretically  predicted  for  a  part  of  the  orbit 

The  measured  acute  polarization  angle  for  the  4  April  i960  lunar 
observation  remained  between  approximately  and  degrees,  as  illustrated  in 
Figure  l8  It  would  appear  from  the  theoretical  rotational  data  that  the  correct 
solution  to  the  ambiguity  problem  is  most  likely  the  0  =  k  -  A0  curve. 

In  order  to  check  the  validity  of  this  assumption,  the  electron 
content  obtained  from  the  various  values  of  0  are  plotted  in  Figure  19>  together 
with  the  theoretical  estimate  Both  the  0  =  r.  +  A  0  and  0  =  2  it  -  A0 
curves  would  indicate  that  the  total  ionospheric  electron  content  was  increasing 
during  the  time  the  electron  contend  below  the  height  of  the  F-layer  maximum  was 
commencing  to  decrease  This  would  require  that  the  electron  production  above 
the  F-layer  maximum  was  continuing  at  a  greater  rate  than  the  electron  loss  below 
the  F-layer  maximum.  The  mechanism  governing  this  physical  process  is  difficult 
to  visualize. 

The  total  electron  content  defined  by  the  0  =  v  -  A  0  curve,  on 
the  other  hand,  closely  resembles  the  variation  of  the  integrated  electron  density 
below  the  F-layer  maximum  It  would  appear  from  the  theoretical  data  that  the 
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Figure  l6.  Average  Faraday  Rotation  from  Trinidad  Radar -Lunar  Observation,  19  March  i960 


RI645 


KTi 

UJ 


Figure  17 .  Comparison  of  the  Measured  Total  ELectron  Content  in  a  One  Square  Centimeter 

Vertical  Column  Through  the  Ionospnere  with  Theoretical  Estimate,  19  March  I96C 
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Figure  ifl.  Average  Faraday  Rotation  from  Trinidad  Radar-Lunar  Observation,  '4  April  1?6G 


iron  Content  in  a  One  Square  Centimeter  Vertical 
eoretical  Estimate,  h  April  I960 


scale  height  for  the  region  above  the  peaK  of  the  F-layer  should  be  somewhat  less 
than  the  80  km  value  employed  in  the  calculations 

The  angular  rotation  data  for  the  complete  lunar  orbit  of  26  Kay  I960 
is  contained  in  Figure  20.  It  is  quite  obvious  that  the  most  acceptable  diurnal 
variation  of  the  polarization  angle  is  given  by  the  0  =  it  -  A  0  ,  B  and  A' 

curves.  There  is  also  excellent  agreement  between  the  experimental  and  theoreti¬ 
cal  toted  electron  contents  as  shown  in  Figure  21 

The  polarization  data  acquired  during  the  whole  lunar  orbit  of  20 
July  i960,  are  given  in  Figure  22  1+  is  seen  "'hat  the  0=  A0  ,  B  and  B’ 

curves  best  fit  the,  theoretical  calculations 

An  examination  of  Figure  23  reveals  that;  during  the  nighttime,  the 
experimentally  determined  electron  content  is  greater  than  the  theoretical  esti¬ 
mates.  As  noted  previously,  this  would  infer  a  scale  height  in  excess  of  80 
kilometers  for  the  region  above  ^he  F-layer  maximum. 

A  composite  of  all  the  total  electron  contents  in  a  cne-square  centi¬ 
meter  vertical  column  through  the  ionosphere  inferred  from  the  lunar  observations 
is  depicted  in  Figure  2k,  There  is  evidence  of  a  diurnal  variation  with  the 

maximum  ionization  content  appearing  in  the  vicinity  of  mid-afternoon:  The  value 
12  2 

of  7  ^  10  electrons/cm^  in  th°  early  hours  of  the  morning  for  the  southerly 

lunar  orbit  of  19  March  i960  is  a  factor  of  about  one-third  times  lower  than  the 

kl 

magnitude  deduced  for  the  northerly  orbits  According  to  Wright  ,  based  on  his 
study  of  electron  density  profiles  in  the  F-layer  along  the  75°w  geographic 
meridian,  the  electron  density  should  maximize  at  about  5°N  geographic  latitude 
which  is  the  region  south  of  Trinidad  Insufficient  radar-lunar  data  recorded 
during  the  time  when  the  moon’s  coordinate  was  at  the  maximum  southerly  declina¬ 
tion  prevent  a  full  analysis  of  this  discrepancy 

From  the  excellent  agreement  found  between  the  electron  content 
determined  from  radar-lunar  observations  and  that  predicted  on  the  basis  of  a 
simple  Chapman  model  for  the  region  above  the  F-layer  maximum,  it  appears  that 
the  ionosphere  may  be  represented,  to  a  first  approximation,  by  a  Chapman  model 
having  a  constant  scale  heignt  cf  80  kilometers  or  less  for  the  daylight  hours 
and  perhaps  80  kilometers  or  more  for  the  nighttime  Wright  has  suggested  that  a 
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Ficure  22.  Average  Faraday  Rotation  from  Trinidad  Radar-Lunar  Observation,  20  July  1060 
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Figure  23.  Comparison  of  the  Measured  Total  Electron  Content  in  a  One  Square  Centimeter  Vertical 
Column  Through  the  Ionosphere  with  Theoretical  Estimate,  20  July  i960 
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scale  height  gradient  of  about  0  2  km/km  for  the  region  above  the  peak  of  the 

F-layer  may  be  more  representative  of  an  ionospheric  model  than  one  with  a 

45 

constant  scale  height  Rocket  sounding  measurement  of  electron  densities, 

20 

as  reported  by  Berning  ,  has  disclosed  a  linear  scale  height  gradient  of  0.22 
km/km  for  the  neutral  gas 


3  Ratio  of  the  Electron  Content  Above  the  F-Layer  Maximum  to  That  Below 


The  diurnal  variation  of  the  ratio  of  the  integrated  electron  density 
above  the  peak  of  the  F-layer  to  that  below,  based  on  all  the  radar-lunar  mea¬ 
surements  recorded  at  Trinidad,  is  summarized  in  Figure  25  It  is  seen  that  the 
ratio,  K,  as  defined  by  Equation  (12),  is  a  variable  function  of  time,  attaining 
a  minimum  of  about  1 . 5  (average  value)  at  approximately  1300  to  1^00  hours  local 
time.  At  night  there  is  a  wide  scatter  of  values  ranging  between  four  and  11. 

It  is  estimated  that  these  results  could  have  an  rms  error  as  high  as  15  per  cent. 


The  ratios  of  the  electron  content  as  measured  by  Bauer  and  Daniels^"" 
compare  rather  favorably  with  the  data  presented  in  this  paper  They  found  the 
ratio  to  be  on  the  order  of  four  to  five  during  three  nights  in  June  before  sun¬ 
rise  and  about  equal  to  three  after  sunrise  In  addition,  for  two  days  in 
November,  the  ratio  was  about  three  both  before  and  after  sunrise. 

4q 

Evans  and  Taylor  have  reported  that,  for  the  months  of  January  and 
February  i960,  the  daytime  ratio,  as  determined  in  England,  appears  to  be  about 
three  while  the  nighttime  ratios  vermin  the  -range. .between- 2  5  sud  10.  Their 
daytime  results  are  slightly  higher  than  those  obtained  at  Trinidad,  but  this 

perhaps  may  be  attributed  to  the  difference  in  geographic  latitude  of  the  obser- 

14  4q 

vation  sites.  E/an’s  early  results  are  in  agreement  with  his  recent  work. 

16 

The  measurements  of  Hill  an!  Dyce  indicate  that  the  ratio  of  the 
electron  content  above  to  below  the  F-layer  maximum  during  the  night  is  on  the 
order  of  1.5  and  during  the  day  about  2. 5 •  It  is  obvious  that  the  nighttime 
result  is  considerably  lower  than  the  Trinidad  data,  This  discrepancy  could  be 
brought  about  by  their  assumption  that  the  function,  f(h)  H  Cos  0,  can  be 
considered  a  constant  in  the  evaluation  of  the  electron  content  integral. 


Equivalent  Slab  Thickness  and  Scale  Height  of  the  Ionosphere 
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Figure  25-  Diurnal  Variation  of  the  Ratio  of  the  Electron  Content.  Above  the  F-layer  Maximum  to  the 
Electron  Content  Below  the  F-layer  Maximum  Based  on  Trinidad  Radar-Lunar  Observations 


A  parameter  presently  employed  to  describe  the  ionospheric  medium  is 
the  equivalent  slab  thickness,  t ,  defined  by  Equation  (13)>  which  is  merely  the 
ratio  of  the  total  integrated  electron  content  in  the  ionosphere  to  the  maximum 
electron  density  at  the  peak  of  the  E-layer . 

A  plot  of  the  diurnal  variation  of  the  maximum  electron  density  of  the 
F-layer,  recorded  at  Puerto  Rico  and  Bogota,  Colombia,  during  the  time  of  the 
Trinidad  lunar  observations,  is  given  in  Figure  26.  The  Puerto  Rico  data  corresponds 
to  the  times  when  the  moon  was  at  its  maximum  northern  declination,  while  the 
Bogota  data  applies  to  the  moon’s  southerly  orbit  of  19  March  i960. 

A  comparison  of  the  total  electron  content  evaluated  from  the  Faraday 
rotation  data,  shown  in  Figure  2k,  with  the  maximum  electron  density  (from  vertical 
incidence  ionospheric  soundings)  Figure  26,  reveals  that  both  parameters  attain 
a  maximum  value  in  the  afternoon. 

It  is  evident  from  Figure  27  that  there  is  a  noticeable  diurnal  influ¬ 
ence  on  the  magnitude  of  the  equivalent  slab  thickness  of  the  ionosphere  (or  F- 
layer ) . 

The  e  pjivalent  slab  thickness  is  a  minimum  between  200  and  370  kilo¬ 
meters  during  the  daytime  and  attains  a  value  as  high  as  800  kilometers  at  night. 

U6 

The  analysis  of  satellite  doppler  radio  signals  by  Ross  and  An derson 

indicates  a  slight  diurnal  influence  on  the  parameter,  t.  The  noon  values  of  r, 

however,  do  shew  a  pronounced  seasonal  variation  which  can  be  fitted  by  a  sinusoidal 

function  of  mean  value  290  kilometers  at  the  equinoxes  and  of  amplitude  100  kilo- 

31  k6 

meters  which  maximizes  near  the  summer  solstice .  J 

With  regard  to  the  presence  of  seasonal  variations  in  the  equivalent 
slab  thickness  of  the  Trinidad  results,  insufficient  data  are  available  to  determine 
this  effect. 

Assuming  that  the  electron  densities  in  the  ionosphere  at  the  time  of 
the  Trinidad  radar-lunar  measurements  were  distributed  in  accordance  with  the 
Chapman  model,  the  scale  height  can  then  be  readily  computed  from  Equation  (14)* 

The  daytime  values  of  scale  height  vary  between  50  and  90  kilometers  with  a  mean 
of  about  75  kilometers,  as  shown  in  Figure  27-  At  night,  there  is  a  wide  scatter 
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Figure  26.  Diurnal  Variation  of  the  Maximum  Electron  Density  of  the  F-layer  During  the  Tine  o 
Trinidad  Radar -Lunar  Observations. 


of  values  ranging  from  50  to  200  kilometers. 

It  should  be  noted  that  the  Trinidad  daytime  evaluations  of  the  scale 
height  of  the  neutral  particles  are  in  close  agreement  with  Kallmann  * s  estimate 
from  rocket  and  satellite  data 

5  Data  Correlation 

An  attempt  was  made  to  determine  whether  there  was  any  correlation 
between  the  calculated  ionosphere  s  parameters,  i  e  ,  the  total  electron  content, 
the  electron  density  ratios,  the  equivalent  slab  thickness  and  the  scale  height, 
and  the  solar  and  geopnysical  conditions  existing  during  the  time  of  the  experi¬ 
mental  measurements 

The  various  solar  and  geophysical  data,  which  were  considered  and  are 
tabulated  in  Table  2,  were  obtained  from  the  bulletins  published  by  the  National 

48 

Bureau  of  Standards 

The  geomagnetic  planetary  index,  k,N,  is  a  mean  three-hourly  magnetic 
reading  based  on  an  arbitrary  scale  of  zero  tc  nine  where  zero  refers  to  a  very 
quiet  day  and  nine  to  an  extremely  disturbed  day  The  relative  sunspot  number  is 
an  index  of  the  activity  of  the  entire  visible  disk  It  is  defined  as 
R  =  K  (lOg  +  s)  where  g  is  the  number  of  sunspot  groups,  s  is  the  total  number  of 
distinct  spots  and  K  is  a  scale  factor  usually  less  than  unity,  its  value  being 
dependent  upon  the  observer  Tne  solar  flux  at  2800  megacycles  is  the  daily 
value  measured  at  noon  in  Ottawa,  Canaia 

Because  of  the  lack  of  sufficient  experimental  data,  this  analysis 
did  not  disclose  evidence  of  any  correlation 

B  CHARACTERISTICS  OF  THE  LUNAR  SURFACE 
1 .  Lunar  Reflection  La ws 

The  supposition  that  the  moon  is  a  rough  body  at  radar  frequencies 
has  led  tc  the  speculation  of  various  scattering  laws  and  functions  to  describe 
the  manner  in  which  radio  waves  would  scatter  from  the  mown  s  surface  A  method 
of  investigating  an  applicable  scattering  law  is  to  characterize  the  decay  rate 
of  the  trailing  edge  of  a  lunar  echo 
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Watts  Meter  (cycles/sec 
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In  oidui  tc  reduce  tht*  e.rfe  l  of  a  no:  i  ieim  fluctuations  m  the 
envelope  of  a  lunar  echo  130  A-scope  photographs,  consisting  of  65  different 
pulse  returns  received  on  the  two  orthogonal  polarizations  were  averaged  to  obtain 
a  representative  sample  for  analysis  purposes  The  resultant  lunar  echo,  shown  in 
Figure  28,  was  obtained  by  dividing  each  pulse  into  C  25-miii i second  intervals. 

At  each  interval  along  rhe  pulse,  an  average  amplitude  was  calculated  from  the  130 
different  data  points 

It  is  seen  that  the  trailing  edge  of  the  pulse  can  be  described  m 
terms  of  the  trigonometric  function.  A  Cot  jj3  id  *■  C  )  J  *  D„  where  the  constants 
A  =  2.5/  B  -  r./55.*  -  ~  '30  and  D  -  -l  i  The  variable  0  is  the  angle  of  incidence 
with  respect  to  the  normal  to  the  surface  of  the  moon  The  relationship  expressing 
0  in  terms  of  the  radio  depth  of  tne  moon  is  contained  in  Append!*  D 


The  power  distribution  of  the  trailing  edge  of  the  average  lunar  echo 
is  compared  with  other  functions  in  Figure  29*  The  third  degree  polynomial, 
a  -  p  9  t  7  O'  ♦  5  0J,  where  a  =  39-06,  p  -  7U.33,  >  =  39-09  and  6  =  -3-38,  also 
accurately  represents  tne  shape  of  the  pulse.  It  is  obvious  however  that  the 

angular  scattering  law,  A,  Gin  2  6/2  0  where  the  constant  A 1  is  equal  to  163, 

II1  u  1 

suggested  by  Leadabrand  does  not  seem  to  apply  to  this  example. 

p 


The  trigonometric  functions,  A  Cos  0  and  Cos  0,  where  A  «  11.2 
and  A~  =  13-65,  also  shown  in  Figure  26,  are  representative  of  the  Lommel -Seel lger 
and  the  Lambert  scattering  lav,  respectively.  When  the  relative  echo  power  is 
plotted  as  a  function  of  the  cosine  of  the  angle  of  incidence,  as  in  Figure  30*  it 
is  found  that  the  Lommel -Seeliger  law  and  the  Lambert  law  reduce  to  straight  lines 
and  that  the  pulse  decay  rate  follows  the  slope  of  the  Lommel -Seel  iger  scattering 
law  displaced  by  approximately  a  factor  of  l/O. 


It  should  be  mentioned  that  Pettengill^  has  reported  that,  following 
the  decay  of  the  initial  specular  component,  the  angular  distribution  of  power 
in  the  moon  echo  obeyed  the  Lambert-type  law  except  for  a  small  amount  of 
limb  brightening  at  the  extreme  ranges.  Since  the  Trinidad  and  Millstone  Hill 
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Figure  30.  Scattering  Laws  Applied  to  Lunar  Echo  as  a  Function  of  the 
Cosine  of  the  Angle  of  Incidence 


investigations  were  conducted  at  approximately  the  same  frequency,  i.e.,  the 
transmission  frequency  at  Millstone  Hill  was  440  megacycles  compared  to  425 
megacycles  at  Trinidad,  the  discrepancy  in  the  results  may  be  attributed  to  the 
fact  that  the  lunar  echo  analyzed  by  Pettengill  was  a  composite  of  24,000  pulses, 
integrated  at  each  500-microsecond  i ange  increment, 

12 

However,  Browne  ,  on  the  other  hand,  has  found  that  an  analysis  of 
lunar  echo  data,  observed  at  a  frequency  of  120  megacycles,  gives  best  agreement 
with  the  Lommel-Seeliger  scattering  model, 

2  Cross  Section  and  Feflertion  Coefficient  of  the  Moon 

In  analyzing  the  cross  sectional  area  of  the  moon  and  the  statistical 
distribution  of  the  amplitudes  of  lunar  signals,  as  discussed  in  the  next  section, 
the  magnitude  of  the  amplitude  of  the  received  echoes  was  measured  at  a  constant 
position  within  the  pulse,  i,e  ,  at  one  millisecond  after  the  beginning  of  the 
pulse  which  corresponds  in  range  to  a  lunar  depth  of  150  kilometers. 

The  cumulative  probability  distributions  of  the  total  cross  section 
of  the  moon  measured  on  12  January  lyQj  during  two  10-minute  periods  are  shown 
in  Figure  31*  The  total  lunar  cross  section  is  merely  the  vectorial  addition  of 
cross  sections  obtained  on  the  orthogonal  polarization  receiver  channels.  Each 
cumulative  distribution  curve  was  calculated  from  over  5000  individual  data 
points.  The  interval  between  l630~l64o  EST  corresponds  to  the  time  when  the 
moon  was  oriented  at  an  azimuth  angle  of  about  72.5  degrees  and  elevation  angles 
of  4.7  -6. 7  degrees,  while  between  2215-2225  EOT  it  was  at  an  azimuth  angle 
varying  from  29  to  15  degrees  and  elevation  angle  of  about  82  degrees. 

It  is  seen  that,  for  the  low  angle  measurements,  50  per  cent  of  the 
total  cross  sections  were  equal  to  or  less  than  8  5  x  lG^  square  meters,  while, 
at  transit,  the  value  reduced  to  6  0  x  10“^  square  meters 

The  lunar  observations  of  2  February  i960,  shown  in  Figure  32,  reveal 
that  50  per  cent  of  the  total  cross  sections  between  l430-l44o  EST  were  equal  to 
or  less  than  5*0  x  lO^1  square  meters,  while,  during  the  interval  2025-2035  EST, 
they  increased  tc  7  5  x  10^  square  meters. 
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Figure  32.  Cunulattve  Probability  Distribution  of  the  Total  Cross  Section  of  the  Moon.  ■  February  i960 


Assuming  that  the  moon  is  a  perfectly  conducting  sphere,  Its  radar 

12 

cross  section  is  then  the  projected  geometric  area  of  the  whole  disk  or  9*5  x  10 
square  meters.  It  follows  from  Equation  (17)  and  Figures  31  and  32  that,  for  50 
per  cent  of  the  observations,  the  power  reflection  coefficient  of  the  moon’s 
surface  appeared  to  lie  between  0.05  and  O.O85.  It  is  estimated  that  the  ex¬ 
perimental  error  incurred  in  this  analysis  should  be  less  than  3  db. 

44 

Measurements  made  by  Frick er  at  a  frequency  of  412  megacycles 
using  a  cw  system,  indicated  a  lunar  power  reflection  coefficient  of  0.07^-  Ac¬ 
cording  to  Blevis  and  Chapman^,  radar  observations  of  the  moon  at  488  megacycles, 
also  made  with  a  cw  system,  revealed  a  power  reflection  coefficient  of  about 
0.05*  The  uncertainty  in  both  experimental  results  are  reported  to  be  less  than 
3  db. 

Since  the  Trinidad  results,  obtained  by  pulsed  radar  techniques  having 
a  pulse  length  less  than  the  radio- depth  of  the  moon,  are  in  excellent  agreement 
with  other  experimental  data  taken  with  cw  systems  in  the  same  frequency  range,  it 
would  appear  that  the  moon  must  be  sufficiently  smooth  for  Fresnel  type  reflec¬ 
tion  to  take  place. 

It  should  be  mentioned  that  values  of  0.1  to  0.15  have  often  been 

quoted  for  the  power  reflection  coefficient  uf  the  moon’s  surface.  The  Trinidad 

4  4  49 

measurements  and  those  reported  by  Fricker  and  Blevis  and  Chapman  would 
indicate  that  these  values  are  somewhat  high. 

3.  Statistical  Analysis  of  Lunar  Echo  Amplitudes 

According  to  theory,  the  amplitudes  of  pulses  reflected  from  an  ir¬ 
regular  surface  which  scatters  electronuignetic  waves  randomly,  i.e.,  with  random 
phases  and  amplitudes,  should  be  Rayleigh  distributed  For  a  surface  which 
contains  both  one  large  smooth  area  and  many  rough  scatters,  the  amplitude  distri¬ 
bution  should  tend  toward  Gaussian 

In  order  to  test  this  hypothesis  for  determining  the  surface  charac¬ 
teristics  of  the  moon,  the  lunar  echo  amplitudes  recorded  in  two  10-rainute  inter¬ 
vals  between  1430-1440  EST  and  2025-2035  EST  on  8  February  i960  were  selected. 

It  should  be  noted  that  the  amplitude  data  from  these  time  periods,  considered 
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in  this  analysis,  are  identical  with  the  lunar  cross  section  data  discussed  in 
Figure  32.  The  statistical  analysis  was  made  on  the  total  signal  amplitude  in 
order  to  reduce  the  effect  of  ionospheric  Faraday  rotation. 

The  probability  density  function  of  the  total  lunar  amplitude,  recorded 
during  moonrise  when  the  libration  fading  rate  was  normally  low,  is  illustrated  in 

Figure  33*  The  theoretical  curves  were  calculated  from  Equation  (27)  with  the 

parameter  defined  by  Equation  (25),  equal  to  28.0,  this  constantbeing  evaluated 
from  the  experimental  data. 

The  presence  of  a  strong  steady  signal  would  be  indicated  by  a 
probability  density  f  motion  which  would  be  of  the  form  of  the  b  =  2  curve.  It  is 

seen  that  the  experimental  points  coincide,  to  some  degree,  with  the  theoretical 

Rayleigh  distribution  (b  o  for  low  amplitude  values,  but  commence  to  diverge 
for  a  normalized  amplitude  of  12  and  greater.  This  is  clearly  evident  if  the 
natural  logarithm  of  the  ratio  of  the  probability  density  function  to  the 
normalized  amplitude  is  plotted  as  a  function  of  the  square  of  the  normalized 
amplitude  as  shown  in  Figure  3^ *  It  is  noted,  from  Equation  (23),  that  for  this 
type  of  a  plot,  the  slope  of  the  Rayleigh  distribution  is  l/2  if  . 

The  experimental  and  theoretical  probability  density  functions  of  the 
total  amplitude  of  the  lunar  echoes  observed  for  10  minutes  near  transit  are 
shown  in  Figure  35-  After  normalization ,  the  parameter,  / ,  for  this  time  period 
was  found  to  be  10.0 

It  is  interesting  to  note  that  the  experimental  data  appear  to  tend 
toward  the  b  ®  1  probability  density  curve  except  for  the  slight  displacement  at 
the  maximum.  It  is  definitely  seen,  from  Figure  36,  that  the  experimental  points 
do  not  correlate  with  the  Rayleigh  distribution  curve. 

A  comparison  of  the  experimental  points  in  Figures  33  and  35  indicate 

kh 

that,  at  high  elevation  angles  where  libration  fading  is  a  maximum,  the  ampli¬ 
tude  signals  were  observed  more  often  having  greater  magnitudes  than  those 
detected  near  moonrise.  This  is  based  on  the  fact  that  the  number  of  pulse 
returns  considered  in  both  analyses  was  the  same. 

An  amplitude  record  of  one  minute  in  duration  (2026:20  -  2027:20  EST) 
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Figure  34.  The  Natural  Logarithm  of  t  .e  Probability  Density  Function  of  the  Total 
Amplitude  of  the  Lunar  Echo,  8  February  i960,  1430-1440  EST. 


PROBABILITY  DENSITY 


31526 


Figure  35. 


The  Probability  Density  Function  of  the  Total  Amplitude 
of  the  Lunar  Echo,  8  February  i960,  2025-2035  EST. 
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Figure  36  •  T5ie  Natural  Logarithm  of  the  Probability  Density  Function  of  the 

Total  Amplitude  of  the  Lunar  Echo,  8  February  i960,  2025-2035  EST. 


taker  irom  the  ^arr-  statist!  '  i.  p  .pulaTicr  was  aiSo  analyzed  to  determine  the 
effect  of  a  shorten*!  $anpl°  on  *bc  experimental  distribution.  In  general, 

it  was  found  tnat  the  amplitudes  of  the  one  minute  sample  vsre  distributed  in  a 
somewhat  similar  fagbton  to  those  of  the  -.0-minute  sample 

The  results  cf  this  analysis  reveal  evidence  of  the  presence  of  a 
slight  specular  component  1  ei  i°:tpd  fr  r  the  i.unar  surface  during  the  time  when 
the  moon  was  near  transiT  This  was  not  tnfr  casr  ,  however,  for  observations 
conducted  at  moon r lee 

I"  possible  that,  instead  i  *:n  large  smooth  region  on  the  sur¬ 
face  of  the  moon,  there  ai  many  smooth  areas*  which  give  rise  to  specular  reflec¬ 
tion  Thai,  this  vojaG  nave  the  eftect  ci  imposing  a  random  pnase  and  random 
amplitude  fluctuation  cn  th:.  rests  i  fan;  ipecular  cumpone^t  It  would  follow  that 
the  steady  signal  which  would  be  noniany  expected  if  the  moon  consisted  of  only 
one  large  smooth  reflection  area,  would,  in  essence,  be  washed  out  or  smeared 

t  Doppler  Fi eqaen^y  Shift 

Because  of  the  relative  motion  of  tbr  moon  with  respect  to  the  earth, 
radio  waves  reflected  from  the  mocn  arc  shined  m  frequency 

The  doppl"-r  frequency  cnifti  wca.  irefi  during  the  partial  lunar  orbits 
of  '  2  January  •+  0  and  I  Tebruarv  1  /£  are  shown  m  Figure  3*7  •  It  is  quite 
evident  that  me  experimental  observations  are  in  agreement  with  the  theoretical 
predict ioHj  Tie  theoretical  curves  w^rr-  computed  according  to  the  method 
proposed  oy  frteker  ,  thr  mathematical  relationships  being  given  by  Equations 
(30  through  (  : 

The  doppler  r.nif:  is  a  maxi  mu m  cf  approximately  +11 30  cps  at  moonrise 
and  a  minimum  at  the  moon  *:  transit  1  -  r.  figure  37 

The  doppler  frequency  shift  at  120  megacycles,  as  reported  by 

12 

Browne  ,  was  on  the  order  of  ±10  cps  fer  observations  taK^n  when  the  hour  angle 
of  the  moon  was  less  than  30  minutes  of  time.  Doppler  shifts  as  large  as  ±1000 
cps  at  >12  megacycles  have  be-^n  recorded  by  Fricber^ 

5  Doppler  Irpquen^ v_3|jpoai 
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Figure  37.  Doppler  Frequency  Shift  of  Lunar  Echoes  Recorded  at  Trinidad 
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Radio  waves  wren  reflected  1  rom  the  moor  experience  a  doppler  spread 
In  frequency  In  addition  to  undergoing  a  doppler  shift  The  former  phenomenon  is 
predominantly  the  result  of  the  moon’s  iibraticn 

uu 

Frizkrr  has  shown  that  the  frequency  spread  at  transit  is  always  a 
maximum,  while  at  moonnse,  it  could  nave  any  value  depending  upon  the  moon’s 
effective  total  li brat Ion  rate  ip  latitude  and  longitude  (sec  Equation  ^hO))- 

7h *  moon  s  Iibraticn  has  1  he  effect  of  imposing  a  rapid  fluctuation 
on  the  amplitude  of  rbe  signal  reflected  from  its  surface  Sample  amplitude  vs 
time  film  records  llluctratmg  the-  presence  of  libration  fading  are  displayed  in 
Figure  36  An  examination  of  Figure  33a  reveals  a  dominant  low  periodic  struc¬ 
ture  which  15  indicative  of  the  minimum  libration  fading  prevalent  at  the  tune  of 
observation  Figure  38c  represents  data  collected  later  in  the  same  day  during 
the  moon's  transit  The  bigr  r  frequency  variation  of  signa)  amplitude,  evident 
in  this  record,  corresponds  to  a  ^ime  of  comparatively  high  libration  rates. 

The  lunar  echoes  obser/ed  at  an  elevation  angle  of  about  16  5  degrees,  Figure 
38b f  show  a  mixture  of  bo^h  low  and  higo  fa ra+es.  It  should  acted  that 
the  relative  difference  of  signal  ampliwudor  in  the  horizontal  and  vertical 
polarizations  from  moonrise  to  -.ran  sit  was  causes.  by  the  change  in  the  rotation 
of  the  plane  or  polarization  by  the  ionosphere.  i  e  ,  Faraday  effect 

In  an  at  temp*  tc  determine  the  correlation  with  tne  ’’effective- 
libration  rate”  doppler,  the  power  spectrum  for  the  two  lunar  observations  were 
obtained  and  compared  with  the  theoretical  doppler  spread 


As  discussed  in  Appendix  E,  the  power  density  spectrum  is  the  Fourier 
transform  of  The  autocorrelation  function  which  is  computed  from  an  amplitude  vs. 
time  function  In  this  analysis,  tve  autocorrelation  function  is  obtained  from 
the  total  lunar  cross  section  which  is  proportional  to  the  square  of  the  signal 
amplitude  The  power  density  spectrum,  discussed  in  this  section,  in  essence 
defines  rhe  frequency -power  content  of  the  time  variation  of  the  total  echoing 
area  of  the  moon 


A  per* ico  of  thr  sample  of  the  8  February  i960  lunar  data  selected 
for  this  analysis  is  shown  ir.  Figure  38a  The  total  signal  power  for  each  pulse, 
i.e  ,  sum  of  the  power  recei  /ed  on  the  two  orthogonal  polarizations,  was 
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ELEVATION  =46.5 


b.  INTERMEDIATE  L  I  BRAT  ION  EFFECTS 
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o. MINIMUM  LI  BRATION  EFFECTS 


Figure  38.  Amplitude  vs.  Time  Film  Records  Shoving  the  Influence 
of  Lunar  Libration,  8  February  i960 
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measured  one  millisecond  ui*er  the  r  ise  f  the  pub'e  Tn is  corresponds  to  the 
coverage  of  a  spherical  cap  on  the  moon  with  a  fractional  radius  of  0  404  (sec 
Appendix  D) 


The  autocorrelation  function  of  the  total  cross  section  for  a  one- 
ninute  interval  at  moon  rise  is  presented  in  Figure  39  The  power  density  spec¬ 
trum  corresponding  to  this  function  is  given  in  Figure  40.  It  is  seen  that  the* 
principle  frequencies  are  all  less  than  one  cps,  i.e  ,  0.27;  0-53;  O.69,  and 
0.9^  with  the  dominant  frequency  being  0.27  eps 


It  should  be  mentioned  mat  a  graphical  Fourier  analyst 
film  record  for  the  same  period  revealed  frequencies  of  0.23;  0-82 > 


51  ^  4Uc 

^  i.  OiiL 

and  1-73  epe 


The  theoretical  maximum  doppler  sp.  ,ud  for  this  particular  time  was 
found  to  be  2  21  cps  for  the  fractional  radius  of  0  404.  This  calculation  wan 
based  cn  the  lunar  libration  rate  constants  presented  in  Table  3  and  applied  to 
Equations  (40)  and  f-ij 

A  similar  analysis  of  the  4  April  \9b0  lunar  data,  also  at  moonriso, 
revealed  the  autocorrelation  function  arid  the  power  density  spectrum  given  in 
Figures  4l  and  42,  respectively  The  principle  frequencies  present  in  the  power 
density  spectrum  are  0,14,  0  27;  0.49;  0  >?,  1  38,  1.62,  1-77;  2-00;  2«30,  2-50; 
2  65,  3  20  and  3  5^  cps.  the  dominant  beirt  0.14  cps  . 

Utilizing  the  libration  nr  t. ant.  tabulated  in  Table  3>  if  Cron  be 
shown  that  the  theoretical  maximum  doppler  spread  for  the  fractional  radius  of 
0  4C4  was  0-66  :ps 


It  is  evident  that,  for  the  two  lunar  observations  that  were 
analyzed,  the  theoretical  maximum  doppler  spreads  exceed  the  dominant  frequency 
of  the  power  density  spectrums  deduced  from  the  experimental  data. 

At  a  frequency  of  423  megacycles,  the  first  Fresnel  zone  on  the  noon 
has  a  fractional  radius  of  approximately  0.007*  This  value  differs  from  the 
fractional  radius  of  0.4o4  employed  in  the  calculations  by  a  factor  of  0.017* 

The  experimental  dominant  frequency  of  the  8  February  i960  observation;  0  27 
is  0  12  times  smaller  then  the  theoretical  estimate  of  2 .21  cps  In  addition. 
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Figrire  k2.  Power  Density  Spectrum  of  Total  Cross  Section,  Moon  Data,  h  April  i960 


the  dominant  frequency  of  O.ih  cps  measured-  on  -  April  I960,  differs  from  the 
theoretical  caicuiar  Oi,  of  0.66  cps  by  a  factor  of  It  appears,  therefor* 

that  while  the  effective  radian  responsible  fci  mod  of  the  amplitude  variation 
is  smaller  than  713  kilometers,  the  radius  of  the  spherical  cap  corresponding  to 
a  fractional  radius  of  OACs  it  is  larger  than  the  radius  of  the  first  Fresnel 
zone,  II/7  kilometers,  by  approximately  10  time. 

As  shown  m  the  photograph  of  the  Irani  portion  of  the  moon  s  face: 
Figure  kf,  there  is  a  comparatively  smooth  area  in  the  center  of  the  moon,  Sinus 
Meaii  large  enough  to  encompass  at  least  10  Fro  no ;  zones.  This  region  is  sur¬ 
rounded  by  a  number  of  mountain  peaks  that  rise  5000  to  3000  feex.^  It  could  be 
possible  that  the  librationai  motion  of  these  peaks  account  for  the  higher  fre¬ 
quency  terms  obtained  in  the  power  density  spectrum.  For  example  the  mountain 
range,  Rhae^icus.  with  a  fractional  radius  of  approximately  0.08-0. 09  could  have 
imparted  at  its  location,  a  fading  rate  of  about  0. 8^-0. 95  ops  on  8  February  i960. 
It  is  interesting  to  note  that  the  p  .  ng  rate,  possibly  due  to  Fhaeticus, 
compares  favorably  with  one  of  tie  frequencies,  1.9^  cps.  present  in  the  power 
density  spectrum  if  that  date.  Thi:  rulanon  is  based  on  Equation  (uOj  for 
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VI.  CONCLUSIONS 


Radar  reflections  from  the  moon  afford  a  method  for  investigating  the  char¬ 
acteristics  of  the  ionosphere. 

It  is  shown  that  the  polarization  ambiguity  problem  of  radar-lunai  reflec¬ 
tions,  prevailing  when  only  one  frequency  is  employed  on  transmission,  can  be 
resolved  by  the  comparison  of  the  experimental  measurements  of  the  angular  rota¬ 
tion  (of  the  plane  of  polarization),  introduced  by  the  ionospheric  Faraday  ef¬ 
fect,  with  theoretical  estimates 

It  appears  that  the  electron  density  distribution  above  the  peak  of  the  F- 
layer  can  be  represented,  to  a  first  approximation,  by  a  Chapman  model  with  a 
constant  scale  height  of  less  than  about  80  kilometers  for  the  neutral  particles, 
during  the  daytime  and  more  than  about  80  kilo.rrrtr^rs  during  the  nighttime  This 
is  based  on  the  nigb  degree  of  correlation  between  the  total  ionospheric  electron 
content  deduced  from  radar-lunar  observations  and  that  predicted  from  ionosonde 
data  and  the  theoretical  mod^l 

The  ratio  of  the  integrated  electron  density  above  F-layer  maximum  to  that 
below  undergoes  a  diurnal  variation,  attaining  a  minimum  average  value  of  ap¬ 
proximately  1.5  between  1300-1400  hours  local  time.  At  night,  there  is  a  scatter 
of  values  ranging  from  four  to  11 

Evidence  is  presented  which  indicates  the  presence  of  diurnal  influence  on 
the  equivalent  slab  thickness  of  the  ionosphere  which  is  defined  as  the  ratio  of 
the  total  ionospheric  electron  content  to  the  electron  density  at  F-layer 
maximum.  The  equivalent  slab  thickness  is  a  minimum  of  200  to  370  kilometers 
during  the  daytime,  while,  during  the  nighttime,  its  magnitude  increases,  attain¬ 
ing  a  value  as  high  as  800  kilometers 

The  scale  height  o^  the  F-region  deduced  from  the  total  electron  content  by 

kl 

the  method  proposed  by  Wright  also  reveals  a  diurnal  influence  The  daytime 
values  vary  between  50  and  90  kilometers  with  a  mean  of  75  kilometers  while 
during  the  nighttime,  the  scatter  is  from  50  to  200  kilometers. 

Lack  of  sufficient  experimental  radar-lunar  data  has  prevented  a  thorough 
analysis  to  determine  the  correlation  of  solar  and  geophysical  phenomena  with  the 
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APPENDIX  A 

IONOSPHERIC  FARADAY  EFFECT 

A.  INTRODUCTION 

When  a  linearly  Cxar‘7ed  electromagnetic  wave  enters  an  ionized  medium  in 
the  presence  of  an  external  magnetic  fi^ld,  such  as  enc our Lered  in  the  ionosjjhere, 
the  wave  separates  into  rvc  independent  components  both,  in  the  general,  case, 
elliptically  polarized  with  opposite  sense  of  rotation.  For  frequencies  in  the 
UHF  range,  the  two  coop  orients  which  are  referred  to  as  the  ordinary  and  the 
extraordinary  waves,  a,e  r_iculariy  .  oi&rizeu. 

As  the  ionosphere  is  traversed,  the  tvc  waves  progress  with  different 
velocities  propagation  which  re si It  in  the  jhase  relationship  between  them  to 
be  continuously  changing.  Or.  leaving  the  lonofphere,  the  circularly  polarized 
components  recombine  t:  form  a  linearly  solarized  wave  which  is  rotated  with 
respect  to  the  original  linear  wave  Thir  phenomenon  is  commonly  referred  to  as 
the  Faraday  effect. 

B .  THEORETICAL  CONSIDERATIONS 

Since  a  circularly  polarized  wave  is  basically  composed  of  two  harmonic 
orthogonal  signals  which  are  ’/•  degrees  cut  of  nase  with  one  another,  only  one 
of  the  linearly  polarized  vftYer  from  each  of  the  "vo  circularly  polarized  components, 
i.e.,  the  ordinary  and  rhe  *  ■/  tr aordinary  waves,  need  be  considered  in  analyzing 
the  polarization  rotational  effects  of  the  ionosphere. 

The  differential  phase  shift,  that  each  of  the  components  undergoes  in  an 
element  of  path  length,  ds  is  given  by 

dp  --  - ds  (A  1) 

O  .A 

and 

dpe  -  — -  d:-  (A-2) 

e 

where  -he  subscripts,  •.  and  .  refer  to  the  ordinary  and  extraordinary  modes  of 
propagation,  respectively,  and  >  and  are  the  wavelengths  associated  with 
the  phase  velocity,  '  and  V  ,  of  the  waves. 

r  TV 
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As  illustrated  in  Figure  A-l,  the  angular  rotation  of  the  plane  of  polarization 
of  the  wave,  &0,  as  it  travels  the  distance,  ds,  is  therefore 


d0  =  1/2  (d0  -  d0  ) 

o  e 


(A-3) 


By  definition, 


Po,e 


o  ,e 


and 


Po,e  n 


(A-U) 


(A-5 ) 


o.e 


where  c  is  the  free  space  velocity;  f,  the  transmission  frequency;  and  n  ,  the 

o ,  e 

refractive  indices  of  the  medium.  It  follows  that  Equation  (A-3)  car.  be  written 
in  the  form 

d0  =  —  An  ds  ( A-6  ) 

where  u  is  the  angular  transmission  frequency  and  An  -  n  -  n  . 

o  e 

This  relationship  defines  the  angular  rotation  of  the  wave  for  a  one-way  propa¬ 
gation  path.  Thus,  the  total  polarization  shift  for  a  two-way  propagation  path 
becomes 

s 


0(s )  =  —  /  An  ds 


(A-7) 


where  s,  and  s  are  the  limits  of  the  path. 
1  2 


In  terms  of  a  vertical  height  variable,  this  expression  is  modified  to 

dh  (A-8) 


h2 


0(h)  =  ~  J  An  f(h) 


9U 


where 


?(h)  - 


r  +  h 

o 


(r  *  h)2  -  (r  Cog  E)‘ 


1/2 


(A-9) 


and  r  Is  the  radius  of  the  earth  and  E  1:  the  elevation  angle  of  the  antenna 
o 

beam. 


According  to  the  magneto- ionic  theory,  l.e  ,  the  theory  of  propagation  of 
electromagnetic  waves  through  an  ionized  medium  an  the  presence  of  an  external 
magnetic  field,  the  complex  index  of  refraction  of  the  medium,  M,  is  given  by  ^  J 


2(  a  -  jP)  - 


(O' 


where 

a 


p  = 


U)  t/ 

U>2 


i  -a  *  03 


:og  e 


(7j) 


(x  -a  ’  J3  ) 


H 


“2  *  U(V 


2 


c 

H  e 
m  c 

4jr  Ne 


An  Q 


1/2 
( A-.10 ) 


( A-ll ) 


ni 


and  wnere  e  is  the  electron  charge  (A. 8  x  10  '  esu):  m ,  the  electron  mass 

pq 

(9*1  x  10  ’  gm);  K,  electron  density  (electrons/cmJ);  v  the  electron  angular 
collision  frequency:  H,  the  earth's  magnetic  field  strength  (gauss);  and  the 
propagation  angle,  i.e.  .  the  angle  between  the  magnetic  field  vector  and  the 
direction  of  propagation.  The  quantity,  He/mc,  is  the  gyromagnetic  frequency  of 
the  electron  about  the  earth’s  magnetic  field  while  Is  the  critical  frequency. 
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or  pla&ma  fr- of  the  '.oniz*  d  ^  U  jn. 

It  is  noted  that  there  are  two  possible  values  for  the  complex  index  of  infrac¬ 
tion  which  indicate  two  different  modes  of  propagation.  The  waves  travel  independently 
In  the  ionized  medium  and  each  has  a  polarization  vector  associated  with  it-  In  the 
general  case,  at  the  low  frequencies,  the  vuves  are  elliptically  polarized  with 

opposite  senses  of  rotation  At  frequencies  on  the  order  of  100  megacycles,  the 

3A 

waves  are  circularly  polarized. 

It  can  be  shown  that,  under  ordinary  ionospheric  conditions,  for  frequencies  in 

2  2  2  2  38 

the  100-mc  range,  it  is  valid  to  assume  J  >  >  v  and  uj  >  w  .  Ihus,  the  complex 
index  of  refraction  reduces  10  a  real  quantity  given  by 


2 


n 

cue 


1  - 


*  H  c  2 

— *  - 7j  Sin  0 

0)  2  uj 

c  c 


uxi)  r*  cj  3  in  0 
H 


2u> 


2 

Cos  6 


1/2 


fA-12) 


where  the  positive  sign  is  associated  with  the  ordinary  mode  of  propagation  while 
the  negative  sign  with  the  extraordinary  mode  of  propagation, 

Assuming  quasi -longitudinal  mode  of  propagation  which  is  represented  by  the 
condition 

2 

4  g  >  Sin  0  Tan  0  (A-13) 

OJ 

H 

it  follows  that  the  difference  In  the  index  of  refraction  between  the  ordinary  and 
extraordinary  waves  is 

An  **  — ~ — r  Cos  6  (A-lU) 

2r:  m  c  f3 

This  approximation  is  valid  at  a  frequency  of  UOQ  m  aacycles  for  values  of  Q  between 

0  <  S  <  89.6 : 


07 


According  to  Equation  (A-8),  the  total  two-way  polarization  rotation  (in 
radians  )  evaluates  to 


0(h  ) 


4.7233  x  10U 


f 


2 


f(h)  H  Cos  0  N  dh 


(A-15) 


where  the  altitude  limits,  and  hg,  are  in  centimeter  units  and  f  in  cps. 

It  is  evident  that,  for  quasi -longitudinal  propagation,  the  magnitude  of 
the  angular  rotation  is  inversely  proportional  tc  the  square  of  frequency  and 
directly  proportional  to  the  magnetic  field  strength. 

The  condition  for  quasi -t ransver se  propagation  is  denoted  by  the  inequality 


2  2 

Sin  0  Tan  Q  >  > 


4'jJ 


U), 


TT 


H 


(a-i6) 


By  simple  manipulation,  it  can  be  shown  from  Equation  (A-12  )  that  the  difference 
Hi  the  refractive  indices  for  the  two  modes  of  propagation  is> therefore,  given  by 


Ue 


a  2 

u 


L  i\ 


3  3  2  - 

6n  m  e  f 


2 

Sm  9 


(A-17) 


Referring  to  Equation  (A-8),  it  follows  that  the  two-way  polarization 
rotation  becomes 


JO 


o-6067  x  10^ 


f(h)  H2 


sin  0 


N  dh 


(A-18) 


It  is  noted  that,  for  the  quasi -transverse  case,  the  angular  rotation  is 
directly  proportional  to  the  magnetic  field  squared  and  Inversely  proportional 
to  the  cube  of  frequency 
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APPENDIX  B 


DETERMINATION  0^  SIGNAL  POWER  REFLECTED  FROM  A  ROUGH  SPHERE 


Any  calculation  oT  the  theoretical  power  reflected  from  a  rough  sphere  must 
include  such  factors  as  the  increment  of  surface  contributing  energy  at  a 
given  instant  of  time  and  a  scattering  la 4  describing  the  surface  characteristics. 

In  order  to  formulate  a  relationship  incorporating  all  the  necessary 
parame ters ,  the  geometric  relation  of  the  incident  pulse  with  respect  to  the 
sphere  must  be  considered. 


A  plane  wave  is  assumed  incident  on  the  surface  of  the  snbere  with  zero 
range  (or  time)  taken  at  the  point  of  tangency.  For  a  pulse  of  xength ,  p  Lin 
units  of  distance  ),  less  than  fhe  radius  of  the  sphere,  r,  "he  incremental  zone, 
olA,  around  the  sphere  :ont ributing  *0  the  reflected  pulse  at  a  given  zone 
depth,  d,  is  as  shown  m  Figure  defined  by  the  angular  limits 


0, 

and 

e_ 


Cos 


i  l  z  (d~r  \  "I 

L  r  J 


co,-L  r 


(  B-l  ) 


1  B-2 ) 


Since  p  -  cl/2  and  d  ct/2, 


U  r 


It  f ol lows,  that 


and 


a 


2 


r 


(B-3) 


(b-M 


where  c  is  the  free  space  velocity;  T,  the  pulce  width  in  unite  of  time;  and  t , 
the  time  associated  with  the  zone  depth. 

For  t  <  T  the  contributions  to  the  received  pulse  originate  within  the 
surface  area  bounded  by  0,-0  and  0?  defined  by  Equation  (B-M 
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RI5II  LEADING  EDGE 

,  OF  PULSE 


D  =  PULSE  WIDTH  (UNITS  OF  DISTANCE) 
d  =  ZONE  DEPTH 
r  =  RADIUS  OF  SPHERE 
dA  =  INCREMENTAL  AREA 

Figure  B-l.  Geometry  of  a  Pulse  Incident  on  a  Rough  Sphere 
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The  method  for  determining  the  magnitude  of  the  signal  reflected  from  a 
rough  sphere>  as  described  bclcw,  is  similar  to  that  developed  by  Grieg,  Metzger, 
and  ,/uer . 

If  P  is  the  power  per  unit  area  incident  on  an  incremental  area,  dA,  of  a 
o 

rough  surface  then  the  effective  incident  power,  ,  is 

PT  P  Cos  0  dA  (B-5  ) 

I  o 

where  0  is  the  angle  between  the  normal  to  the  surface  and  the  direction  of 
the  incident  wave.  The  angle,  0^  ,  in  Figure  B-l  corresponds  to  0  in  this 
expression . 

If  the  surface  is  not  perfectly  reflecting,  the  reflected  power,  p  _s  given 

R 

by 


=  o  p  Cos  0  dA  ( B-6  ) 

R  o 

where  p  is  the  power  reflection  coefficient  of  the  surface. 

It  is  convenient  to  define  P,T  &s  the  power  per  unit  solid  angle  reflected 

n 

normal  to  the  area  dA  For  Lambert  scattering,  the  power  radiated  m  any 

direction  is  proportional  to  the  cosine  of  the  angle  between  the  scattered  ray 

and  the  normal  to  the  surface.  Thus,  the  power  scattered  per  unit  solid  angle 

is  simply  P  Cos  7  where  y  is  any  angl^  measured  from  the  normal  to  dA. 

N 

The  unit  sol^d  angle  enclosed  by  dA,  as  shown  in  Figure  B-2,  is  given 

by 


r 


It  can  be  readily  shewn  fron;  spherical  geometry  that 

dA  =  2ji  r  3in  7  d>  IB-8; 
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i r.  follow  "-iv  fhe  total  power  intercepted  (or  reflected)  by  the  upper 
hemi$phe“<'  of  Figure  E-2  i:  therefore, 
r./2 


P„  =  2n  P 


II 


r 

o 


Shi  7  Co:  7  cLy 


(B-9) 


or  simply 


P  «  It  P 
p  :j 


(B-10) 


.'■Vt  .  rf  i  let  tec  i;  wei  i  1  defined  ly  Equation  (B-6)  :  t  is  evident 


that 


pn  = 


p  P  Cos  0  dA 
o 


(B-ll ) 


Taking  into  account  Laxuber"  fr  law  once  again  it  follows  that  the  power 
per  ’ir.it  colli  angle  s*  uttered  back  in  the  direction  of  0,  i  e  ,  incident 
direction  becomes 


p  P  Cos  6  dA 
o 


( B-12 ) 


It  is  the  transmitted  power  and  G  is  the  gain  of  the  transmitting  antenna 
relative  to  an  isotronic  antenna,  then 


P  Cn 


P 

o 


t-n  p 


(B-13) 


where  P  is  the  distance  from  the  transmitting  source  to  the  sphere* 


Ire  pr- intercepted  at  the  receiving  antenna,  p  ,  is  slirply 


A 

pr  ‘  pe  (  —%—) 

R 


(B-14) 


where  A  is  the  effective  antenna  area  defined  by 


(B-15) 


•r 


where  G  is  the  gain  of  the  receiving  antenna  and  s  Is  the  transmitted  wavelength 
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Since  the  element  of  area,  dA,  can  be  expressed  in  terms  of  the  angle  9  by 
(assuming  01  =  Q  in  Figure  (B-i)  ) 

dA  =  2n  r2  Sin  9  &Q 

on  substituting  Equations  (B-12),  (B-13)*  (B-15),  and  (B-l6)  in  Equation 

it  is  seen  that  the  total  received  power  for  a  monostatic  system  can  be  expressed 

by  the  integral 


P  = 
r 


2  £ 

Pt  0  *  °8 


y  r- 


S 

8 j  Sin  0  Cos  0  10 

0. 


where  is  the  echoing  area  of  a  smooth  reflecting  body  defined  by 


2 

a  -nr  p 
s  H 


( B- 18  } 


It  is  noted  that  the  terms  in  the  parenthesis  of  Equations  (b-i7*  describe 
the  general  radar  equation,  i  e.,  the  power  intercepted  fro  wThic 

reflects  radiation  isotropically. 

On  integrating  between  the  limits  j  =  0  and  0o  =  Jt/2,  it  is  found  that 
the  total  power  received  from  a  rough  sphere,  scattering  in  accordance  with 
Lambert's  law,  is  8/3  times  greater  than  that  from  a  smooth  body 


With  regard  to  a  rough  surface  obeying  the  Lotmel- See lager  scattering  law 

the  power  scattered  per  unit  solid  angle,  P  ,  is  independent  of  the  scattering 

f» 

angle,  y ,  i.e.,  the  angle  between  the  scattered  ray  and  the  normal  to  the 
surface.  Thus,  the  total  power  intercepts d  by  the  upper  hemisphere,  shew,  in 
Figure  B-2,  becomes 

*/2 

r 

PR  -  2 it  P^  j  Sin  7  dy  (B-19  } 

o 


or 


PR  *  2*  P 


N 


( B-20 ) 


1C  4 


It  can  be  readily  shovm  that  the  power  received  from  a  rough  sphere 
evaluates  to 


P  = 


2  2 

Pt  G  A  °S 


(4rt  )"  R 


k 


Si  a  0  Cos  0  ^0 


(B-21) 


Integration  between  the  limits  G1  ~  0  and  0 0  -  n/2  reveals  that  the  power 

scattered  from  a  Lommel-Seeligcr  type  sidierical  surface  is  twice  that  of  a 
smooth  sphere 

The  theoretical  shape  of  pulses  reflected  from  a  rough  moon  are  depicted  in 
Figure  B-3*  The  relative  power  scale  has  been  normalized  with  respect  to  the 
Lommel-Seeliger  lav  It  is  seen  that,  for  short  pulse  lengths,  more  energy 
is  scattered  back  from  the  limb  of  the  moon  for  the  Lornmel-SeeJ  lger  scattering 
case  than  for  Lambert  scattering  On  tui.>  pui  ticular  plot,  a  smooth  moon  would 
reflect  a  constant  relative  power  of  0.5>  irrespective  of  pulse  width. 


Figure  B-3.  Shape  of  pulses  Reflected  From  a  Rough  Toon 


APPENDIX  C 


ORBIT  OF  A  CELESTIAL  BODY 

The  relationship  of  the  celestial  coordinate  system  with  the  geograpnic 
coordinate  system  _s  illustrated  in  Figures  C-l  and  C-2.  The  declination,  Df 
of  a  point  on  the  celestial  sphere  is  analogous  to  geographic  latitude  since  it 
is  the  angular  displacement  from  the  equatorial  plane  measured  positive 
northward.  The  right  ascension,  RA,  whi^h  corresponds  to  geographic  longitude 
is  the  angular  position  about  the  celestial  pole  It  is  measured  positive  east 
from  the  meridian  containing  the  Vernal  Equinox  (or  the  First  Point  of  Aries) 
which  is  the  intersection  of  the  plane  of  the  ecliptic,  i.e.,  the  plane  of  the 
sun’s  orbit  about  t^e  earth,  with  the  celestial  equator.  The  Greenwich  Hour- 
Angle  of  Aries,  GHA^,  is  the  angular  distance  measured  positive  westward  along 
the  celestial  equator  from  the  Greenwich  meridian  to  the  Vernal  Equinox,  The  local 
hour  angle,  0,  is  the  angular  displacement  between  the  meridians  containing  the 
observation  site  and  the  point  of  the  celestial  sphere,  and  is  measured 
positive  westward. 

It  can  be  readily  shown  by  applying  the  law  of  cosines  to  the  spherical 
triangle  of  Figure  C-3  that  the  function  relating  the  geographic  and  celestial 
coordinates  is  given  by. 

Sin  D  =  Sin  E  Sin  L  +  Cos  E  Cos  L  Cos  A  (C-l) 

where  L  is  the  latitude  of  the  observation  sight,  and  A  and  E  are  the  azimuth 
and  elevation  angles  of  the  celestial  body  as  observed  at  the  site,  respectively 
For  the  special  case  when  A  =  l80°,  this  equation  simplifies  to 


E  =  — —  +  D  -  L  (C-2 ) 

For  A  -  0° 

E  =  L  ±  (  — | -  D)  (C-3) 


where  the  two  solutions  exist  only  when  the  condition,  L  >  |^(rt/2) 
satisfied. 


-  D 
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D  - 
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DECLIN  ATI  Of 
RIGHT  ASCE 
POINT  OF  V 


oordinate  System 


VERNAL  EQUINOX 
(ARIFS) 


Figure  ’-2. 


;eles  ial  Time  Coordinates 


Figure  C-3 •  Astronomical  Configuration 
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According  to  Equation  (Ol),  it  is  seen  that  the  azimuth  and  elevation 
orientations  of  a  celestial  body  are  only  dependent  upon  the  observer’s  latitude 
and  the  deciinatior *  Since  the  declination  of  a  star  is  relatively  constant, 
its  orbit  as  observed  from  a  fixed  location  on  the  earth's  surface,  remains  the 
same  from  day-to-day  and  is  symmetrical  about  the  observer’s  meridian.  With 
regard  to  the  sun  and  the  moon,  however,  both  orbital  paths  are  not  constant 
since ‘their  declinations  are  variable  functions  of  time. 

The  duration  of  time  ^hat  a  celestial  body  will  remain  above  the  horizon  at 
any  latitude  can  be  readily  determined.  Referring  to  Figure  C-3>  it  can  be 
shown  from  spherical  trigonometry  that,  in  terms  of  units  of  time,  the  azimuth  and 
elevation  coordinates  can  also  be  expressed  by 

Sin  E  ~  Cos  0  Cos  L  Cos  D  ^  Sin  L  Sin  D  (C-^) 

and 

Tan  A  = - Sl“  - -  (C-5) 

Tan  D  Cos  L  -  Sin  L  Cos  0 

I 

It  should  be  noted  that  the  local  hour  angle,  0,  is  expressed  in  sidereal  time 
units  and  not  in  solar  time.  Since  the  sidereal  day  is  shorter  than  the  solar 
day,  a  star  will  cross  an  observer’s  meridian  approximately  three  minutes,  55*91 
seconds  earlier  every  day  . 

The  time  at  which  a  celestial  body  such  as  a  star  crosses  a  given  angular 
position  can  be  determined  by  utilizing  Equations  (C-l)  and  (C-5  )  In  conjunction 
with 


0  =  3o0’  *  C.HA*  -  RA  -  X  (C-6) 

where  >  is  the  longitude  of  the  observer  measured  positive  westward.  This 

y 

expression  is  derivable  from  Figure  C-2.  Tne  parameter,  GHA  ,  is  listed  as  a 
function  of  time  of  day  in  the  Nautical  Almanac . 

With  regard  to  determining  the  time  of  passage  of  the  moon  or  the  sun, 

Equation  ;--6)  is  modified  to  the  form 

0  *  Gill  -  >  (C-7) 

^  w 
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where  CHA,  the  Greenwich  Hour  Angle,  is,  in  essence,  the  local  hour  angle  Uus 
angle,  which  is  also  tabulated  in  the  Nautical  Almanac,  is  measured  positive 
westward  from  Greenwhich  meridian. 
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APPENDIX  D 

THE  GEOMETRY  OF  THE  MOON 


According  to  the  theory  of  Fresnel  diffraction  as  pertaining  to  physical 

56 

optics  ,  the  radius  of  the  n*tb  Fresnel  zone,  a  is  expressed  by 


where  X  is  the  wavelength  of  the  transmitted  monochromatic  light  and  S  and  P 
are  the  distances  between  the  diffracting  aperture  and  the  light  source  and  the 
observing  screen^  respectively 

For  a  source  at  infinite  distance,  the  dimension  of  the  n‘th  zone  reduces  to 


a  * 
n 


nip 


1/2 


f  D-2  ) 


With  regard  to  the  concept  of  specular  reflection  from  a  smooth  spherical 
surface  as  applied  to  the  radar  problem,  the  radius  of  the  ncth  zone  is  modified 
to 


n  X  R 

a  ■=  — - - 

n  L  2 

This  relationship  implies  that  the  transmitter  and  receiver  are  situated  at  the 
same  location  and  that  the  target  dimensions  are  large  compared  to  the  wave¬ 
length  of  the  incident  radiation 

A  plot  of  the  radius  of  the  Fresnel  zones  on  the  surface  of  the  moon  at  a 
frequency  of  ^25  megacycles  is  given  in  Figure  D-l.  The  average  earth-moon 
separation,  F,  is  assumed  to  be  385^000  kilometers  for  this  calculation* 

The  fractional  radius  oi  the  moon,  K,  also  depicted  in  Figure  D-l,  is  defined 
by 

J  a 

K  *=  -2—  (D-M 

where  r  is  the  radius  of  the  moon, 

» 


l/i 


D  \ 
J  J 


113 
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For  an  electromagnetic  wave  incident  on  the  moon,  the  depth  of  penetration 
into  the  moon  is  given  by 


where  c  is  the  free  space  velocity  and  t  is  the  two-way  travel  time 

Ac  shown  in  Figure  D-2,  the  radius  of  the  base  of  the  spherical  cap,  a  , 
defined  by  the  penetration  depth,  is  expressible  by 


a  *  2  rd-d 


1/2 


(D-6) 


The  angle  of  incidence,  0,  with  respect  to  the  normal  to  the  lunar  surface 
is  related  to  the  fractional  radius  by 


0 


-1 

=*  Sin 


K 


(D-7) 


A  plot  of  the  angle,  9,  as  a  function  of  the  travel  time  of  the  wave,  the 
penetration  depth  and  the  fractional  radius  of  the  moon  is  given  in  Figure  D-3* 
A  photograph  of  the  face  of  the  moon  illustrating  the  functional  relationship 
between  these  various  parameters  is  shown  in  Figure  D-U. 
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APPENDIX  E 

AUTOCORRELATION  FUNCTION  AKD  POWER  DENSITY  SPECTRUM 

A  statistical  or  random  function  of  time,  y(t),  such  as  random  noise 
current,  is  one  in  which  the  value  of  the  function,  y(t),  does  not  depend  in  a 
definite  manner  on  the  independent  variable  time.  The  function  can  best  be 
described  in  terms  of  a  statistical  or  probability  distribution. 

The  autocorrelation  function  is  a  measure  of  the  average  correlation 
between  two  values  of  the  same  time  series.  It  is  similar  to  the  correlation 
coefficient,  often  used  in  statistical  analysis,  which  is  a  measure  of  the 
linearity  between  two  variables. 

The  autocorrelation  function  is  defined  by 

/T/2 

y(t )  y  (t  +  t )  dt  (E-i ) 

-  T/2 


For  .a  random  process  which  is  stationary  in  time,  i.e.,  the  causes  of  the 
ran  iom  fluctuations  do  not  change  with  time  so  that  the  process  is  independent 
of  the  time  of  origin,  \jr(t)  is  an  even  function  with  respect  to  t* 


The  autocorrelation  function,  in  normalized  form,  can  also  be  written  as 


■  y(t)  -  yjt) 

*(t)  -  - - —— 

y(t)  -  yTTT 


Some  of  the  properties  of  the  normalized  autocorrelation  function  are  namely; 

t(C)  a  1;  'I'(t)  £  i>(o)  ;  and  for  a  random  phenomenon  without 

the  presence  of  a  steady  signal,  Jt(t) — »  0  as  t-><» 

The  power  density  spectrum,  u)(f)  of  the  function  y(t)  is  the  Fourier  transform 
of  the  autocorrelation  function  In  other  words  , 
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and  similarly  the  autocorrelation  1  re  ion  is  uei  .v  from  the  Courier 
transform  of  the  power  density  spectrum. 


In  simpler  fora;  the  relationship  between  J/(t)  s--d  oj(f  .an  be  expressed 


i(-t)  -•  2  J  -j(f)  Cos  2nf-r  df 


(E-4) 


and 


Ji£) 


r 


J 


f (t  )  Cos  2n ft  It 


With  reference  to  the  analysis  if  the  amplitude  versus  time  functior.  of 

lunar  echoes,  the  resultant  signal  reflected  from  the  surface  of  the  moon 

undergoes  a  rapid  fluctuation  due  to  libration*  Since  the  lunar  libration  rate 

changes  continuously  throughout  the  orbit,  a  statistical  description  of  the  _  . , 

/  v 

amplitude  record  obtained  at  one  particular  time,  is  therefore,  only  valid  for  1 
one  portion  of  the  orbit  „ 


Radar  data  are  usually  acquired  at  equally  spaced  time  intervals  which  are 


dependent  upon  the  pulse  repetition  frequency  of  the  transmitted  signal  The 
autocorrelation  function  for  this  fcm  c  ’  a  discrete  time  series  can  be  computed^y^ 


by 
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n-rh 


*<r>"  ~ —  z  xi  xi+„ 

n-rh  L — 1 


1-0 


X  is  the  time  scries,  and  r  **  0.1 
n 

interval  A  «  h  At  for  an  integer  h>o  where  Ar 


where  X  ,  X  , 
o’  i 


m  where  nh 
t '  T*e  interval 


adjacent  values. 


The  computation  of  the  power  density  spe  trum  can  be  repi  :>:■  '•  \  by 

m-1 

tfo)  +  2  V  +(i  )  cos.  •— 

L _ i 

irxl 


V(r )  - 


+  ^m)  Cos  rrr 


where  V(r)  is  the  raw  estimate  of  the  saoothi.  nower  density  at  the  nc  ual 
frequency  r/(2  m  At) 


9  > 


